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ABSTRACT 


A  computer  analysis  has  been  developed  to  predict  the  behavior 
of  tall  buildings  subjected  to  lateral  loads.  The  buildings  considered 
are  composed  of  shear  walls  and  frames.  Each  shear  wall -frame  unit 
undergoes  equal  lateral  deformations  and  the  structure  as  a  whole  does 
not  deform  torsionally.  The  analysis  is  performed  on  a  subassemblage 
consisting  of  a  lumped  shear  wall  system  and  a  simple  substitute  frame. 
The  load-deflection  relationship  for  the  structure  is  determined  up  to 
the  ultimate  load.  The  formation  of  plastic  hinges  in  the  substitute 
frame  has  been  considered  as  well  as  the  inelastic  action  of  the  wall. 

In  addition,  the  secondary  moments  produced  by  the  axial  forces  in  the 
columns  acting  through  the  sway  displacements  of  the  structure  are 
considered . 

The  lumping  procedure  used  to  form  the  subassemblage  is 
discussed  and  the  analysis  of  the  subassemblage  in  both  the  elastic  and 
inelastic  range,  including  the  P  -  A  effect,  is  presented.  The 
results  of  the  analysis  are  compared  with  those  obtained  by  other 
methods.  The  behavior  of  a  shear  wall-frame  structure  is  discussed  and 
illustrated  by  an  example. 
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CHAPTER  I 


INTRODUCTION 

Many  types  of  buildings  are  composed  of  shear  walls  coupled  to 
structural  frames.  The  shear  walls  may  be  on  the  exterior  face  of  the 
building  or  in  the  form  of  interior  cores  for  elevators  or  stairs. 

Shear  walls  generally  have  stiffnesses  which  may  be  much  greater  than 
those  of  the  members  composing  the  structural  frame.  It  has  been  common 
to  design  the  frame  to  resist  the  vertical  loads  acting  on  the 
structure  and  the  walls  to  resist  the  lateral  forces  The  strength 

and  stiffness  of  the  shear  wall  was  estimated  by  the  same  methods  used 
for  frame  members. 

The  above  procedure  is  acceptable  for  low  buildings,  but  as 

more  slender  structures  are  considered,  the  interaction  between  the 

shear  wall  and  the  frame  becomes  significant.  If  the  frame  alone 

resisted  the  lateral  forces  it  would  deflect  as  shown  in  FIGURE  1.1a, 

in  this  case  the  sway  deformations  of  each  story  are  approximately  equal. 

On  the  other  hand,  the  free  shear  wall  would  deflect  as  shown  in 

FIGURE  1.1b.  Since  the  two  systems  are  connected  the  redistribution  of 

forces  will  result  in  a  final  deflected  shape  as  shown  in  FIGURE  1.2. 

The  frame  system  will  restrain  the  wall  in  the  top  stories  and  as  a 

result  forces  will  be  induced  in  the  frame  that  cannot  be  predicted  by 

(2) 

assuming  the  total  lateral  load  to  be  resisted  by  the  wall 


1 


2 


Since  the  structure  under  consideration  is  highly  indeter¬ 
minate,  methods  which  account  for  the  inelastic  action  will  more  reliably 
predict  the  behavior  of  the  system.  In  tall  buildings,  the  axial  loads 
in  the  columns  are  significant  and  produce  additional  overturning  moments 
commonly  referred  to  as  P  -  A  moments;  where  P  is  the  axial  load  and 
A  is  the  displacement  due  to  sidesway.  The  effect  of  these  "secondary" 

moments  is  to  reduce  the  capacity  of  the  structure  to  resist  lateral  loads. 

(3) 

These  moments  must  be  considered  in  the  analysis  of  tall  building  frames 
The  use  of  computers  has  increased  considerably  during  the  past 
decade,  and  this  makes  it  possible  to  analyze  complicated  problems. 

Computer  methods  are  particularly  suitable  for  elastic-plastic  analysis, 
since  the  structure  must  be  analyzed  in  several  stages  of  deterioration, 
in  order  to  obtain  the  complete  response. 

The  above  discussion  clearly  indicates  the  importance  and 
complexity  of  combined  shear  wall -frame  systems.  The  present  investi¬ 
gation  attempts  to  predict  the  load-deflection  behavior  of  these  systems, 
both  in  the  elastic  and  the  inelastic  range.  The  P  -  A  moment  is 
included  in  the  analysis  and  the  corresponding  reduction  in  the  lateral 
load-carrying  capacity  of  the  structure  is  examined.  Example  structures 
were  analyzed  and  the  results  are  compared  to  those  obtained  by  other 


methods . 
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FIGURE  1.1  FREE  DEFLECTED  SHAPE  OF  FRAME  AND  SHEAR  WALL  FIGURE  1.2  FINAL  SHAPE  OF  COMBINED  SHEAR 


CHAPTER  II 


PREVIOUS  INVESTIGATIONS 

Recently  the  behavior  of  structures  composed  of  shear  walls 

(4  5  6  7  8  9) 

and  rigid  frames  has  aroused  a  great  deal  of  interest  v  *****  '. 
Methods  of  elastic  analysis  have  been  presented  but  no  real  attempt  has 
been  made  to  study  the  behavior  of  the  system  in  the  inelast  .  range. 

The  traditional  methods  of  frame  analysis  make  no  provision 
for  wall  type  elements  having  a  finite  width.  The  girders  attached  to 
the  walls  will  be  subjected  to  moments  due  to  the  vertical  displacement 
of  the  connection  points.  These  moments  are  in  addition  to  those 
which  would  be  computed  by  assuming  the  wall  to  be  a  line  element. 

Also,  the  equations  resulting  from  traditional  analyses  of  structures 
composed  of  stiff  elements  connected  to  flexible  members  are  ill- 
conditioned  ^0)  so  that  the  resulting  solutions  may  be  inaccurate. 
However,  such  structures  have  been  analyzed,  and  no  difficulties 
reported  ^  ^ . 

Coupled  shear  walls  may  be  analyzed  by  replacing  the  discrete 

(12  13  14) 

connecting  beam  by  a  continuous  lamina  of  equivalent  stiffness  v  ’  ’  ' 

However,  in  order  to  obtain  a  closed  form  solution,  it  is  necessary  to 
assume  constant  structural  properties  throughout  the  height  of  the 
building.  An  iterative  solution  has  been  developed  for  this 
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problem,  which  is  capable  of  accounting  for  changes  in  configuration  of 

(4)  (21 

the  shear  walls  .  This  procedure  is  similar  to  that  developed  v  '  for 

the  solution  of  combined  shear  walls  and  frames  and  will  be  discussed  in 

detail  below. 

Expressions  for  the  lateral  stiffness  of  various  types  of 

( 5l 

beam  and  column  arrangements  have  been  developed  by  Cardan  v  .  Each 
story  is  isolated  by  assuming  that  the  column  has  an  inflection  point 
at  mid-height.  The  moments  and  forces  developed  by  the  various 
structural  arrangements  are  then  assumed  to  be  continuously  distributed 
over  the  height  of  the  shear  wall.  A  differential  equation,  expressing 
the  equilibrium  of  the  shear  wall  is  formulated  by  assuming  the 
structure  to  have  constant  properties  throughout  its  height.  This  last 
assumption  makes  the  method  unsuitable  for  most  practical  problems. 

Parme  ^  has  combined  the  resistance  of  the  shear  wall  with 
that  of  the  frame  to  obtain  an  expression  for  the  total  lateral  load  on 
a  floor  as  a  function  of  the  stiffness  of  members  (including  the  shear 
wall)  and  the  displacements  and  rotations  of  the  floor  considered  and 
the  adjacent  floors.  The  response  of  the  shear  wall  is  formulated  in 
terms  of  finite  difference  equations.  Finally,  the  solution  of  the 
resulting  set  of  simultaneous  equations  determines  the  lateral  displace¬ 
ment  of  each  floor.  The  analysis  neglects  the  axial  deformation  of  the 
columns  and  the  effect  of  the  wall  width  on  the  moments  in  the  girders 
adjacent  to  the  wal 1 . 
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FIGURE  2.1  shows  the  model  used  by  Gould  Each  story  is 

represented  by  rigid  bars  connected  by  extensional  and  rotational  springs. 
The  individual  story  properties  and  the  effect  of  adjacent  floors  are 
considered  in  determining  the  lateral  resistance  of  a  particular  story. 

The  analysis  assumes  constant  story  heights  and  constant  shear  wall 
stiffness  and  neglects  the  effect  of  axial  deformations  of  the  columns 
and  shear  deformations  of  the  members.  The  lateral  deformations  of  the 
wall  are  represented  by  finite  difference  expressions.  This  procedure 
leads  to  a  set  of  simultaneous  equations,  one  for  each  story. 

As  described  above,  the  matrix  methods  have  also  been  used 
by  Clough,  King  and  Wilson  to  analyze  frame  shear  wall  structures. 
These  procedures  include  the  effects  of  wall  width  and  axial  distortions. 
No  attempt  was  made  to  perform  an  inelastic  analysis  or  include  the 

secondary  moments  due  to  the  P  -  A  effect. 

(8) 

Bandel  v  '  expressed  the  deformed  shape  of  the  combined  truss 
(or  wall)  and  frame  system  by  using  a  power  series.  Expressions  for  the 
total  internal  and  external  work  of  the  system  are  derived  and  the 
constants  of  the  actual  deformation  functions  are  determined  from  the 
solution  of  a  set  of  simultaneous  equations.  The  equations  result  from 
the  minimization  of  the  total  potential  energy  of  the  system.  It  is 
necessary  in  this  case  to  replace  the  shear  wall  by  an  equivalent  truss 
system. 

(q) 

Rosenbleuth  and  Holtz  v  '  assumed  that  an  initial  configuration 
may  be  obtained  from  the  differential  equation  solution,  by  assuming  the 
structure  to  be  completely  uniform.  This  initial  configuration  is  then 
improved  in  later  cycles  of  iteration  by  considering  the  actual  properties 


of  the  structure. 
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f  2) 

Khan  and  Sbarounis  v  '  have  developed  an  iterative  method,  in 
which  the  building  is  lumped  to  form  a  structural  model  (FIGURE  2.2). 

The  model  consists  of  a  shear  wall  and  a  simplified  frame  system.  The 
elastic  analysis  of  the  present  work  is  based  on  a  similar  approach. 

The  method  accounts  for  the  variation  in  wall  and  frame  properties. 

The  effects  of  base  rotation,  plastic  rotation  of  the  wall  and  secondary 
deflections  in  the  frame  have  also  been  discussed.  Influence  curves  are 
provided  to  estimate  the  distribution  of  shear  between  the  wall  and  the 
frame  for  a  wide  range  of  structural  proportions. 

The  above  discussion  summarizes  the  existing  literature  dealing 
with  the  analysis  of  shear  wall-rigid  frame  structures  in  the  elastic 
range.  So  far  no  attempt  has  been  made  to  determine  the  inelastic 
behavior  of  such  structures.  Literature  is  available,  however,  concerning 
the  inelastic  behavior  of  unbraced  multi -story  planar  frameworks  acting 
without  shear  walls. 

Several  theories  are  available  to  predict  the  behavior  of 
multi -story  frames.  In  order  to  correlate  the  results  of  the  analysis 
with  test  results,  however,  certain  factors  must  be  accounted  for.  As 
the  load  is  increased  on  the  frame,  portions  of  the  member  will  yield, 
thus  decreasing  the  stiffness  of  the  frame  under  additional  loading. 

For  tall  slender  structures  the  secondary  moments,  produced  by  the  axial 
column  loads  acting  through  the  sway  displacements,  reduce  the  capacity 
of  the  frame  to  resist  lateral  loads.  In  addition,  for  these  structures. 
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the  redistribution  of  moment  caused  by  column  shortening  may  become 
significant.  It  has  been  found  that  an  elastic-plastic  analysis  which 
includes  the  above  effects  will  closely  predict  the  response  of  the 
structure  up  to  the  ultimate  load. 

The  second-order  elastic-plastic  analysis  of  the  complete 

(15  16  17) 

structure  can  be  performed  using  the  displacement  method  '  *  *  . 

A  second-order  elastic  analysis  is  used  up  to  the  formation  of  the  first 
hinge.  The  stiffness  is  then  reduced  by  inserting  an  idealized  hinge  in 
the  structure  at  this  point.  The  analysis  is  then  continued  for  succes¬ 
sive  increments  of  load.  At  each  load  increment,  this  procedure  results 
in  a  set  of  simultaneous  equations  which  are  the  equilibrium  equations 
for  the  deformed  structure.  The  joint  rotations,  vertical  deflections, 
and  story  translations  are  determined  from  the  solution  of  the  above 
set  of  equations.  The  end  forces  in  the  members  are  then  determined  by 
substitution.  Each  time  a  new  hinge  is  formed,  the  determinant  of  the 
stiffness  matrix,  before  and  after  hinge  formation,  is  calculated.  When 
the  determinant  passes  from  a  positive  value  to  zero  or  a  negative  value, 
the  maximum  capacity  of  the  structure  has  been  attained. 

Isolated  parts  of  a  multi -story  frame  have  been  analyzed  as 
(3) 

independent  units  v  .  The  points  of  inflection  in  the  columns  are 
assumed  to  be  at  mid-height,  so  that  two  imaginary  horizontal  slices 
through  the  structure  isolate  one  floor  level  with  its  attached  columns 
as  shown  in  FIGURE  2.3.  Subsequently  this  model  is  reduced  to  four 
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subassemblages,  each  consisting  of  a  column  restrained  by  the  adjoining 
beams.  The  load-deflection  curve  for  each  subassemblage  is  obtained  by 
performing  an  inelastic  analysis;  the  load-deflection  relationship  for 
a  particular  story  is  then  obtained  by  summing  the  responses  of  the 
individual  subassemblages. 

In  the  case  of  frame-shear  wall  structures,  the  behavior  of  a 
portion  of  the  shear  wall  cannot  be  isolated.  The  deflections  and 
rotations  of  each  part  of  the  wall  depend  on  the  deformation  of  the  wall 
as  a  whole.  Thus,  a  subassemblage,  which  would  include  the  influence  of 
the  shear  wall,  cannot  be  obtained  in  a  rational  manner. 
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FIGURE  2.1  GOULD'S  MODEL 
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FIGURE  2.2  KHAN'S  MODEL 


FIGURE  2.3  CENTRAL  PORTION  OF  TWO  STORY  UNIT 
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ANALYTICAL  MODEL  FOR  THE 

FRAME-SHEAR  WALL  STRUCTURE 

In  order  to  reduce  the  analysis  of  a  complete  structure  to 

manageable  terms,  the  actual  structure  is  replaced  by  a  simplified 

(18) 

structural  model.  Lightfoot  v  '  has  used  a  similar  scheme  to  replace  a 
multi-story  frame,  such  as  that  shown  in  FIGURE  3.1a,  by  the  equivalent 
frame  shown  in  FIGURE  3.1b.  The  substitute  frame  consists  of  columns 
having  stiffnesses  equal  to  the  sum  of  the  stiffness  of  all  the  actual 
columns  in  a  particular  story,  restrained  at  each  floor  by  a  beam, 
having  a  stiffness  equivalent  to  the  sum  of  the  stiffnesses  of  all  the 
actual  beams  at  that  floor  level. 

In  FIGURE  3.1a  and  3.1b,  K-^,  ,  Kgg  etc.  are  the  beam  end 

moments  developed  for  a  unit  rotation  at  each  end  of  the  beam.  For 
prismatic  members, 
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In  the  above  expression  I  and  L  represent  the  moment  of 
inertia  and  length  respectively  of  beam  12. 

The  equivalent  beam  restraint,  K^E  ,  at  each  floor  is 

given  by 


2  (K12  +  K 
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+ 


k34) 


and  the  stiffness  of  the  equivalent  column,  ,  is  given  by 


KCE  K1  +  K2  +  K3  +  K4 

where  K-j ,  K^,  and  are  the  column  stiffnesses  in  the  story. 

The  equivalent  frame  can  then  be  analyzed  under  the  given 
external  loads.  The  moments  in  the  actual  frame  can  be  determined 
from  the  analysis  of  the  substitute  frame,  since  the  moments  and 
forces  in  each  member  are  proportional  to  the  ratio  of  the  stiffness 
of  the  actual  member  to  the  stiffness  of  the  equivalent  member. 

(19) 

Frischmann  and  Prabhu  v  '  have  adopted  a  similar  method  to 
lump  the  member  stiffnesses  of  multi -story  framed  structures  and  have 
also  used  the  method  to  analyze  inter-connected  shear  wall  structures. 

In  the  present  analysis,  the  actual  shear  wall-rigid  frame 
structure  has  been  lumped  into  a  structural  model  consisting  of  two 
systems:  the  Frame  System  and  the  Shear  Wall  System.  These  are 
shown  in  FIGURE  3.2.  In  FIGURE  3.2,  EL  F.  is  the  equivalent  beam 
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representing  the  actual  beams  which  span  between  two  columns; 

B.  W..  is  the  equivalent  beam  representing  the  actual  beams  which 
link  the  columns  to  the  shear  wall  at  the  i  th  floor  of  the 
structure.  In  the  model  the  shear  walls  in  the  actual  structure  have 
been  replaced  by  a  single  shear  wall  with  an  equivalent  stiffness  and 
width.  Similarly  the  columns  in  the  actual  structure  have  been  lumped 
to  form  an  equivalent  column.  In  the  analysis  of  the  model,  joints 
B..  and  F.  are  subjected  to  equal  rotations  and  F.  is  free  to 
translate  laterally.  The  rotational  base  springs  for  the  column  and 
the  shear  wall  reflect  the  foundation  conditions. 

Building  frames  are  made  up  of  various  combinations  of 
frames  and  shear  walls.  FIGURE  3.3  shows  the  plan  of  a  building 
consisting  of  two  exterior  frames  and  a  central  shear  wall.  The  lateral 
loads  are  assumed  to  be  applied  in  a  direction  parallel  to  the  exterior 
frames.  Assuming  the  floors  to  be  infinitely  rigid  in  their  own  plane, 
so  that  all  points  on  a  floor  are  forced  to  undergo  equal  lateral 
translations,  the  building  in  FIGURE  3.3  can  be  replaced  by  the  system 
shown  in  FIGURE  3.4.  It  has  been  assumed  that  the  floor  slabs  offer  no 
rotational  restraint  to  the  shear  wall,  so  that  the  system  resisting 
the  lateral  loads  consists  only  of  the  two  exterior  frames  linked  to  the 
shear  wall.  In  FIGURE  3.4,  Kg  and  Kg  are  the  stiffness  (EI/L)  of 
each  individual  column  and  beam  in  the  building. 


- 
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The  system  in  FIGURE  3.4  can  be  replaced  by  the  equivalent 
system  shown  in  FIGURE  3.5.  In  any  story,  the  columns  are  lumped  into 
a  single  equivalent  column.  The  frame  system  consisting  of  the 
equivalent  column  restrained  by  a  beam  is  linked  to  the  wall  system. 

In  FIGURE  3.5,  it  has  been  assumed  that  both  ends  of  the  left  beam  will 
undergo  equal  rotations.  Since  the  restraint  offered  by  the  beam  must 
be  equivalent  to  the  restraints  on  the  column  shown  in  FIGURE  3.4,  the 
stiffness  of  the  left  beam  in  FIGURE  3.5  is  given  a  stiffness  equal  to 
twice  the  sum  of  the  stiffnesses  of  all  the  beams  in  the  actual  structure. 

In  other  types  of  buildings,  beams  may  directly  connect 
adjacent  columns  to  the  shear  wall  system.  Generally  these  beams  offer 
relatively  little  restraint  to  the  lateral  deflection  of  the  structure 
and  a  simple  and  conservative  way  of  forming  the  analytical  model  would 
be  to  neglect  the  action  of  these  members  and  lump  the  remaining  beams 
in  the  left  beam  as  before.  This  would  neglect  the  restraint  offered 
by  the  beams  or  slabs  connected  directly  to  the  shear  wall.  If  this 
restraint  is  significant,  the  problem  of  finding  an  equivalent  system 
for  such  a  building  becomes  complex  and  a  complete  solution  has  not  been 
attempted. 
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(a)  ACTUAL  FRAME  (b)  EQUIVALENT  FRAME 


FIGURE  3.1  LIGHTFOOT'S  SUBSTITUTE  FRAME 
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FIGURE  3.2  ANALYTICAL  MODEL 
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FIGURE  3.3  PLAN  VIEW  OF  EXAMPLE  BUILDING 
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FIGURE  3.4  LUMPED  EXTERIOR  FRAMES 
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FIGURE  3.5  LUMPED  ANALYTICAL  MODEL 


CHAPTER  IV 


ANALYSIS  OF  MULTI-STORY 

FRAME-SHEAR  WALL  STRUCTURES 

In  the  present  analysis,  it  is  assumed  that  the  behavior  of 
the  shear  wall-rigid  frame  building  can  be  represented  by  the  substitute 
structural  model  described  in  the  previous  chapter.  Each  floor  of  the 
building  is  assumed  to  form  a  diaphragm  of  infinite  rigidity.  Thus  all 
points  on  a  particular  floor  will  have  the  same  lateral  deflection. 
Torsional  effects,  due  to  non-coincidence  of  the  centre  of  load  and  the 
centre  of  resistance  of  the  building  have  not  been  considered. 

It  is  assumed  that  lateral  buckling,  lateral-torsional  buckling 
and  local  buckling  of  the  members  are  prevented;  only  the  in-plane 
behavior  of  the  model  is  considered.  The  effects  of  axial  shortening 
and  shear  deformation  of  the  members  and  the  effect  of  axial  load  on  the 
stiffness  of  the  columns  have  been  neglected.  In  the  columns,  it  is 
assumed  that  idealized  plastic  hinges  will  form  at  the  ends  of  the 
members  only. 

The  moment-curvature  (M-0)  relationships  for  the  individual 
members  of  the  frame  are  assumed  to  be  ideally  elastic-plastic  as  shown 
in  FIGURE  4.1.  In  FIGURE  4.1,  M^c  represents  the  plastic  moment 
capacity  of  the  member,  reduced  for  axial  load  (in  the  columns)  while 
0pC  represents  the  curvature  corresponding  to  MpC  ,  assuming  ideally 
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elastic  behavior. 

The  moment-curvature  relationship  for  the  shear  wall  system  is 
read  into  the  computer  program  as  a  series  of  moment  and  curvature 
co-ordinates.  FIGURE  4.2  shows  a  typical  plot  of  the  moment,  M  , 
non-dimensional ized  as  M/M  ,  versus  the  curvature,  0  ,  non-dimension¬ 
al  i zed  as  0/0^  ;  where  is  the  plastic  moment  capacity  of  the  wall 

and  0^  is  the  elastic  curvature  correspondi ng  to  (0^  =  M  /El) 

For  the  examples  that  follow  only  two  branches  of  the  M-0  curve  have 
been  used. 

The  loading  is  assumed  to  be  applied  statically.  The  structure 

is  analyzed  under  a  constant  vertical  load  and  monotonical ly  increasing 

lateral  load.  The  analysis  begins  by  applying  the  entire  lateral  load 

to  the  shear  wall.  The  corresponding  shape  of  the  deflected  shear  wall 

is  shown  by  the  dashed  line  in  FIGURE  4.3.  The  deflection  at  the  i  th 

level  is  ^  where  i  represents  the  floor  under  consideration 

and  the  superscript  (1)  represents  the  initial  iteration  cycle.  The 

frame  is  forced  into  this  deflected  shape  and  the  corresponding  (negative) 

frame  shears  are  computed.  The  shear  wall  is  reanalyzed  by  applying  the 

negative  shears  from  the  frame  analysis.  The  resulting  deflected  shape 

is  shown  by  the  solid  line  in  FIGURE  4.3.  The  net  deflection  from  the 

(2) 

previous  positon  of  the  shear  wall  is  Aw  .  '  .  The  frame  is  again 

forced  into  the  corresponding  deflected  shape  and  the  negative  shears 

recomputed.  In  the  actual  process  used,  the  frame  is  not  forced  into  the 

(2) 

deflected  shape  characterized  by  Aw  . v  '  but  instead  a  forcing  formula 

is  applied  to  speed  convergence.  The  new  position  is  shown  by  the 

(3) 

broken  line  and  is  characterized  by  the  deflection  A .  '  .  The 
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above  process,  and  the  necessary  modifications,  are  described  in  detail 

in  the  following  paragraphs. 

A  computer  program  has  been  developed  to  assist  in  the 

analysis  of  shear  wall-frame  structures.  The  program  takes  into  account 

yielding  in  the  frame,  the  inelastic  action  of  the  shear  wall  and  the 

secondary  moments  caused  by  the  P-A  effect.  The  axial  force,  P  ,  on 

the  structure,  acting  through  a  sidesway  displacement,  A  ,  causes  an 

overturning  moment  in  addition  to  that  produced  by  the  lateral  loads, 

known  as  the  P-A  moment.  To  account  for  this  effect,  the  equilibrium 

equations  are  formulated  on  the  deformed  structure.  This  "secondary" 

moment  causes  significant  reductions  in  frame  strength  and  must  be 

(3) 

included  in  the  analysis  of  tall,  slender  structures 

The  method  used  to  analyze  the  model  consists  of  a  series  of 
iterations.  The  steps  involved  are  described  below: 

(A)  The  vertical  load  acting  on  the  frame  is  computed.  This  load 
remains  constant  while  the  lateral  load  on  the  structure  is 
incremented  to  determine  the  load-deflection  relationship. 

The  initial  value  of  the  lateral  load  is  selected  so  that  the 
structure  remains  elastic. 

(B)  Each  story  of  the  shear  wall  is  divided  into  segments.  The 
entire  lateral  load  is  applied  to  the  shear  wall  and  the 
moment  at  the  centre  of  each  segment  is  calculated  by 
considering  the  wall  as  a  free  cantilever.  The  rotations 
and  deflections  are  calculated  at  every  segment  using  moment- 
area  principles.  The  vertical  displacement  at  the  junction 
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of  the  shear  wall  and  the  frame  is  computed  as  the  product  of 

the  rotation  of  the  wall  and  one-half  the  wall  width.  FIGURE 

4.4  shows  the  deflected  shape  of  the  shear  wall  under  the 

entire  lateral  load,  where  0, ,  .  ^  ,  A,,  .  ^  and 

w  i  W  i 

6^  .  ^  are  the  rotation,  lateral  displacement  and  vertical 
displacement  of  the  wall.  The  initial  subscript  W  represents 
the  shear  wall,  while  i  represents  the  floor  under  considera¬ 
tion  and  the  superscript  (1)  represents  the  initial  iteration 
cycle.  The  story  height  is  Hs  . 

(C)  The  frame  system  is  forced  into  the  deflected  shape  of  the  wall 

as  defined  by  the  deformations  computed  above.  This  step  is 
shown  in  FIGURE  4.5. 

Next  the  frame  joint  rotations  are  computed.  In  FIGURE  4.5, 

B.  is  the  beam-to-column  joint,  W.  is  the  joint  at  the 
connection  of  the  wall  and  the  right  side  beam  and  F.  is  the 
left  end  of  the  left  side  beam.  Under  lateral  loads  beam 
B.  F.  will  have  a  point  of  contraf lexure  approximately  at 
midspan.  This  is  enforced  in  the  structural  model  by  assuming 
that  the  joints  B.  and  F.  have  equal  rotations  and  that  the 
left  end  support,  F.  ,  is  free  to  translate. 

At  each  joint  the  slope-deflection  equations  are  written  and 
by  using  the  moment  equilibrium  equation,  an  expression  is 
developed  which  relates  the  joint  rotation  to  the  deformations 
of  the  adjacent  joints  and  the  member  properties.  The  detailed 
derivation  of  this  equation  is  given  in  APPENDIX  A. 


V  ' 
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(D) 


(E) 


A  set  of  equations  is  generated  by  applying  the  moment 
equilibrium  equation  to  all  the  beam-to-column  joints  in  turn. 
The  joint  rotations  are  then  computed  using  the  Gauss-Seidel 
Iteration  Method.  The  moments  at  potential  hinge  locations  in 
the  members  of  the  substitute  frame  can  now  be  computed.  Next 
the  vertical  and  horizontal  shears  at  the  junctions  of  the  wall 
and  frame  system  are  computed. 

At  each  floor  level,  the  moment  to  be  applied  to  the  wall 
system  for  the  next  cycle  of  iteration,  .  ,  is  the 

algebraic  sum  of  the  moment,  M^B  .  ,  at  the  end  W.  of  the 

right  beam  and  the  moment  produced  by  the  beam  shear  multiplied 
by  half  the  wall  width.  FIGURE  4.6  shows  typical  forces  at 
floor  level  i.  In  the  initial  iteration,  the  moment,  Mtl  .  , 

is  assumed  to  be  zero.  This  is  corrected  in  subsequent  cycles. 

The  net  out-of-balance  forces  are  applied  to  the  wall  and  the 

(2)  (2) 
resulting  deflections.  A-.  ,  rotations  0^  .  ' 


and  vertical  displacements  6^ 
in  Step  (B). 


(2) 


are  computed  as  described 


(F) 


In  order  to  speed  convergence,  the  frame  is  not  forced  into  the 

(2) 

deformation  mode  defined  by  .  ,  etc.  Instead  the  frame 


is  forced  into  the  position  defined  by: 


6W  i 


(2)  _  0W  i 


(1) 
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A', 


A  (1)  A  (1) 

(2)  =  AW  i  '  AW  i 

"a  (1)  a  (2) 

flW  i  -  AW  i 


.  (4-2) 


and  6'  .  (2)  =  0'  .  (2)  .  D,.  ./2 

W  l  W  l  W  l 


(4-3) 


where  D 


W  i 


is  the  width  of  the  shear  wall  at  the  i  th 


floor.  The  analysis  now  returns  to  Step  (C)  and  the  steps 
(C)  to  (F)  are  repeated.  For  any  cycle  (after  the  initial 
cycle),  the  forcing  Equations  (4-1),  (4-2),  and  (4-3)  are 
replaced  by: 

„  0)  a  (n-1) 

(n)  _  °W  i  '  ®W  i  . (4-4) 

d  II*  “  - 

1  0  (n-1}  -9  (n) 

0  W  i  0  W  i 

A  (1)  A  (n"1} 

Ai  (n)  =  AW  i  •  AW  i  . (4-5) 

Wi  "a  (n-])-A  (n) 

AW  i  AW  i 

and  ^'w  i  ^  =  01 W  i  ^  ^  ’  DW  i^2  . (4-6) 


where  n  denotes  the  cycle  in  progress. 


For  any  cycle  (except  the  initial  cycle)  convergence  tests 
for  rotations  and  deflections  are  performed  at  the  end  of 
Step  (D)  to  determine  whether  further  iteration  is  required. 


The  system  has  converged  if  the  quantities 
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W  i 


W  i 


and 


are  less  than  a  specified  limit 


A 


W  i 


(normally  0.01). 

If  the  system  has  not  converged  as  defined  above,  then  force 

F.  and  moment  M,ID  .  (FIGURE  4.6)  at  each  floor  are  applied 
i  WB  i  rr 

on  the  shear  wall  and  the  Steps  from  (E)  to  (G)  are  repeated. 

On  the  other  hand,  if  the  system  has  converged,  the  P-A 
effect  is  included.  The  total  axial  load  at  each  floor  level 
of  the  structure  has  been  read  into  the  program.  The  increased 
story  shears  due  to  the  P-A  effect  are  simulated  by  increasing 
the  applied  lateral  loads.  The  additional  force  .  is 
calculated  on  the  basis  of  the  deformations  obtained  from  the 
anlaysis  under  the  lateral  loads  H.  .  The  derivation  is 
given  in  APPENDIX  B.  Next  the  lateral  loads  (H.  +  H-j  . )  are 
applied  on  the  shear  wall  system  and  Steps  (B)  to  (G)  are 
repeated.  The  deformations  of  the  structure  obtained  from  this 
analysis  are  compared  with  the  deformations  obtained  from  the 
analysis  under  lateral  loads  H.  only. 

If  the  deformations  obtained  from  the  two  analyses  do  not  agree 
(within  a  specified  limit)  new  additional  forces  .  are 
computed  as  above  using  the  deformations  of  the  structure  taken 
from  Step  (H).  Steps  (B)  to  (I)  are  then  repeated. 


If  the  deformations  agree  the  shear  resisted  by  the  frame  is 
then  deducted  from  the  total  lateral  load  to  determine  the  shear 
to  be  resisted  by  the  wall. 
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The  formation  of  hinges  in  the  structure  will  change  the  moment 
equilibrium  equations  at  each  beam-to-column  joint.  FIGURE  4.7 
shows  the  possible  hinge  locations  between  stories  i-1  and 
i+1  .  The  modified  moment  equilibrium  equations  are  developed 

in  APPENDIX  A. 

At  the  end  of  each  load  increment  the  deflected  shape  of  the 
structure  is  obtained  as  well  as  the  distribution  of  forces  and 
moments.  In  order  to  determine  the  load-deflection  curve  for 
the  structure,  the  lateral  load  is  incremented  and  the  Steps 
(B)  to  (L)  are  repeated. 

For  each  increment  of  lateral  load,  the  structure  is  treated  as 
a  separate  problem.  The  manner  of  including  the  inelastic 
action  of  the  structure  implies  that  the  stiffness  deterioration 
in  the  analysis  lags,  by  one  load  increment,  the  actual 
stiffness  deterioration.  To  minimize  this  effect,  large  load 
increments  are  used  in  the  elastic  range  and  smaller  increments 
in  the  inelastic  range  of  behavior. 

The  analysis  was  performed  using  a  Fortran  Program  written  for 
the  I.B.M.  360  Computer.  By  slightly  modifying  the  basic  program  it  is 
possible  to  use  the  program  for  elastic  first  order,  elastic  second  order 
or  inelastic  analysis.  The  printout  of  the  Fortran  Program  has  been 
included  as  APPENDIX  C. 

In  the  'MAIN1  program,  the  number  of  problems  to  be  solved  is 
read  in.  The  flow  chart  is  shown  in  FIGURE  4.8.  Then  Subroutine  'SR' 
is  called.  The  flow  chart  for  Subroutine  'SR'  is  shown  in  FIGURE  4.9. 


. 
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This  subroutine  reads  in  the  data  required  for  the  analysis  of  the  first 
structure.  It  computes  the  co-ordinates  of  each  floor  from  the  base  and 
the  stiffness  of  all  the  members.  All  the  input  quantities  are  printed 
out.  Control  is  then  returned  to  the  'MAIN'  program. 

The  maximum  number  of  load  increments  to  be  used  in  the 
analysis  is  specified.  The  initial  load  system  (i.e.  first  increment) 
is  applied  to  the  shear  wall  system. 

Subroutine  ' BAKA '  is  called.  The  flow  chart  is  given  in 
FIGURE  4.10.  'BAKA'  computes  the  distribution  of  lateral  load  between 
the  frame  and  the  shear  wall  systems.  The  moments  and  deformations  at 
each  segment  of  the  wall  are  computed.  The  convergence  formula  is  then 
applied  (except  in  the  first  cycle)  to  the  deflections  and  rotations. 

The  deformations  as  computed  by  the  convergence  formula  are  enforced  on 
the  frame  system.  The  joint  rotations  of  the  frame  system  are  computed 
in  Subroutine  'FRAME'  (FIGURE  4.11)  by  the  Gauss-Seidel  Iteration  Method. 
In  'BAKA',  the  moments  at  all  locations  are  compared  with  the  plastic 
moment  capacities  (except  during  the  first  load  increments).  If  a  hinge 
is  detected,  the  moment  is  set  equal  to  the  plastic  moment  capacity  for 
all  subsequent  load  increments.  The  deformations  are  compared  with 
those  of  the  previous  cycle  of  iteration.  If  the  deformations  are 
within  1  %o  of  one  another,  control  returns  to  the  'MAIN'  program.  If 
the  deformations  are  not  within  the  1  %  limit,  the  process  is  repeated 
within  'BAKA1.  In  this  manner  the  forces  and  deformations  (elastic)  for 
a  given  loading  condition  are  completely  known. 

Equivalent  horizontal  forces  simulating  the  P-A  effect  are 
computed.  These  are  added  to  the  initial  horizontal  load  and  the  total 
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forces  are  then  applied  to  the  structure.  The  structure  is  analyzed  by 
recalling  Subroutine  ' BAKA ' .  The  deflections  obtained  from  the  two 
analyses  (i.e.  with  and  without  axial  load)  are  compared;  if  they  are 
within  1%  ,  the  results  of  the  last  analysis  are  printed  out  by 

Subroutine  ' ROFA'  (FIGURE  4.12)  and  'SR  3'  (FIGURE  4.13). 

On  the  other  hand,  if  the  deflections  are  not  within  1%  , 

new  equivalent  horizontal  loads,  are  computed  on  the  basis  of  the  latest 
set  of  deflections.  The  process  is  repeated  until  convergence  is 
achieved. 

Inelastic  action  of  the  shear  wall  is  considered  in  the  'MAIN' 
program.  At  each  segment  of  the  shear  wall,  the  computed  moment  is 
compared  with  the  plastic  moment  capacity.  If  the  moment  exceeds  the 
plastic  moment  capacity  of  the  story,  the  moment  of  inertia  of  the 
plastified  wall  segment  is  reduced  according  to  the  M-0  diagram.  This 
reduced  stiffness  is  used  in  subsequent  steps  in  the  analysis. 

In  Subroutine  'SR  1'  (FIGURE  4.14),  the  formation  of  plastic 
hinges  in  the  members  of  the  frame  system  is  detected.  If  a  hinge  is 
detected,  the  location  of  the  plastic  hinge  is  printed  out  with  the 
magnitude  of  the  moment.  Control  then  returns  to  'MAIN'  program. 

If  the  structure  is  still  elastic,  the  horizontal  load  is 
incremented  and  the  structure  is  reanalyzed.  On  the  other  hand,  if  a 
hinge  has  formed  at  some  location  in  the  structure,  the  horizontal  load 
is  decremented  so  that  the  structure  is  still  elastic.  The  structure 
is  then  reanalyzed  with  this  reduced  load.  Subsequent  increments  in 
load  are  reduced  in  magnitude  to  trace  the  inelastic  range  as  closely  as 
possible.  The  program  stops  when  the  deformations  of  the  structure  are 
unusually  large.  It  then  returns  to  beginning  to  analyze  the  next  structure. 
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FIGURE  4.1  IDEALIZED  MOMENT-CURVATURE  DIAGRAM  -  FRAME  MEMBERS 


1.0  <*>,/<*>„ 


FIGURE  4.2  IDEALIZED  MOMENT-CURVATURE  DIAGRAM  -  SHEAR  WALL 


WALL 
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FIGURE  4.4  SHEAR  WALL  UNDER  LATERAL  LOAD 
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FIGURE  4.5  DEFLECTED  FRAME 


WALL 
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FIGURE  4.6  FRAME  FORCES  AND  MOMENTS  ON  SHEAR  WALL  FIGURE  4.7  POTENTIAL  PLASTIC  HINGE  LOCATIONS 
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FIGURE  4.8  MAIN  PROGRAM 
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FIGURE  4.8  MAIN  PROGRAM  (CONTINUED) 
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FIGURE  4.9  SUBROUTINE  SR 
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FIGURE  4.10  SUBROUTINE  BAKA 
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FOR  INSUFFICIENT 
NUMBER  OF 
CYCLES 


FIGURE  4.10  SUBROUTINE  BAKA  (CONTINUED) 
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FIGURE  4.11  SUBROUTINE  FRAME 
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FIGURE  4.12  SUBROUTINE  ROFA 


FIGURE  4.13  SUBROUTINE  SR  3 
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SR  1 


FIGURE  4.14  SUBROUTINE  SR  1 


CHAPTER  V 


RESULTS  AND  DISCUSSION 

The  method  of  analysis  described  above  has  been  used  to  analyze 
various  multi-story  structures.  In  this  section,  the  results  of  the 
analyses  of  two  frames  are  discussed.  The  first  example  is  used  to 
demonstrate  the  validity  of  the  analysis  in  the  elastic  range  while  the 
second  illustrates  the  influence  of  several  factors  on  the  behavior  of  a 
typical  multi-story  structure  both  in  the  elastic  and  the  inelastic  range. 

The  first  frame  was  analyzed  primarily  in  order  to  provide  a 
check  on  the  analysis  in  the  elastic  range.  This  frame  had  been  analyzed 
previously  in  REFERENCE  2.  The  actual  building  is  a  fourteen-story 
structure,  rectangular  in  plan,  with  nineteen  bays  of  11 1 -6 "  in  the  long 
direction  and  three  bays  of  20' -0"  in  the  short  direction.  The  building 
has  been  analyzed  for  a  uniformly  distributed  lateral  load  applied 
perpendicular  to  the  long  side  of  the  building. 

The  geometry  of  the  equivalent  structure  is  given  in  FIGURE  5.1. 

The  equivalent  stiffness  of  the  left  beam  is  twice  the  sum  of  the  stiffnesses 
of  the  three  beams  spanning  20'-0"  at  a  given  floor  level.  The  equivalent 
moments  of  inertia  of  the  various  members  (including  the  shear  wall)  have 
been  listed  in  TABLE  5.1.  The  moments  of  inertia  of  the  actual  members 
have  been  given  in  REFERENCE  2.  The  structure  was  analyzed  elastically 
by  a  different  lumping  procedure  in  REFERENCE  2.,  however,  the  model  used 
in  this  report  represents  a  further  simplification  in  that  it  contains  a 
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single  equivalent  columa  The  modulus  of  elasticity  for  both  the  wall 

2 

and  the  frame  was  taken  as  3000  k/in  .  The  rotational  spring 
constants  at  the  base  of  the  column  and  the  shear  wall  were  set  at 
relatively  large  values  so  that  the  fixed  end  condition  is  approached  to 
represent  the  original  structure. 

The  distributed  lateral  load  was  replaced  by  single  concentrated 
loads  acting  at  each  panel  point.  In  order  to  obtain  values  for  the 
vertical  loads  acting  on  the  structure,  reasonable  uniformly  distributed 
live  and  dead  loads  were  assumed.  The  values  of  the  total  vertical  load 
acting  at  each  level  of  the  structure  are  listed  in  TABLE  5.1.  These 
values  are  required  for  the  second-order  elastic  analysis. 

First  order  and  second  order  (including  the  P-A  effect) 
elastic  analyses  were  performed  on  the  structure.  The  results  of  the 
first  order  elastic  analysis  were  compared  with  the  results  presented  in 
REFERENCE  2.  The  agreement  was  satisfactory .  FIGURE  5.2  shows  a  graph 
of  the  story  number,  N  ,  versus  the  proportion  of  the  total  base  shear 
carried  by  the  frame  system,  V^/V^  .  The  solid  lines  indicate  the 
results  of  the  first  order  elastic  analysis  using  the  model  shown  in 
FIGURE  5.1  and  the  dashed  lines  represent  the  results  obtained  from 
REFERENCE  2.  The  ratio  V.p/Vb  and  thus  the  actual  shear  in  the  story,  is 
relatively  constant  except  near  the  bottom  of  the  structure.  Since  the 
applied  shear  increases  from  the  top  down,  the  top  stories  of  the  frame 
are  forced  to  carry  shears  which  are  much  greater  than  the  applied  shears. 
This  is  due  to  the  pull  exerted  on  the  top  stories  of  the  frame  by  the 
stiff  shear  wal 1 . 

The  deflected  shape  of  the  fourteen  story  structure  is  shown  in 
FIGURE  5.3.  This  figure  plots  the  story  number,  N  ,  versus  the  lateral 
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deflection,  A  ,  at  the  various  floor  levels.  The  deflections  plotted 
correspond  to  a  lateral  load  of  20  psf  lumped  at  each  floor  level.  The 
solid  line  represents  the  deflected  shape  obtained  from  a  first  order 
elastic  analysis.  The  dotted. line  shows  the  deflected  shape  which  results 
when  the  P-A  effect  is  taken  into  account.  Although  the  P-A  effect 
is  not  too  significant  in  this  example,  it  becomes  important  for  more 
flexible  structures. 

In  this  example,  the  shears  developed  in  the  top  stories  were 
considerably  higher  than  those  applied.  To  study  the  effect  of  the  wall 
stiffness  on  the  shear  distribution,  six  additional  analyses  were 
performed.  Various  combinations  of  shear  wall  stiffnesses  were  considered 
in  the  top  four  stories  of  the  structure;  the  properties  of  the  lower 
stories  remaining  constant.  In  TABLE  5.2,  the  moments  of  inertia  of  the 
top  four  stories  of  the  shear  wall  are  listed  along  with  the  corresponding 
shears  resisted  by  the  frame  system.  In  this  table  the  frame  shears  are 
listed  as  a  percentage  of  the  applied  story  shears,  V^/V  .  In  all 
cases,  the  percentage  story  shear  resisted  by  the  top  story  of  the  frame 
was  decreased  when  compared  with  the  results  of  the  original  analysis. 

On  the  other  hand,  the  percent  shear  resisted  by  the  lower  stories  of  the 
frame  was  increased.  For  example,  considering  cases  (4)  and  (5)  in 
TABLE  5.2,  there  is  a  decrease  in  in  the  13th  and  14th  story  and 

a  correspondi ng  increase  in  the  11th  and  12th  story. 

The  second  example  considered  is  a  twenty-four  story,  three  bay 
frame.  This  frame  was  designed  using  both  the  allowable  stress  and  the 
plastic  strength  techniques  in  REFERENCE  20.  The  general  layout  of  the 
frame  is  shown  in  FIGURE  5.4a.  The  analytical  model  used  in  the  present 
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case  is  shown  in  FIGURE  5.4b.  The  moments  of  inertia  and  plastic 
moment  capacities  of  the  beams  used  in  the  model  structure  are  listed 
in  TABLE  5.3.  The  total  vertical  loads  acting  on  the  structure  and  the 
properties  of  the  columns  are  listed  in  TABLE  5.4.  The  bracing  system 
shown  in  FIGURE  5.4a  was  replaced  by  a  shear  wall  as  shown  in  FIGURE  5.4b. 
No  attempt  was  made  to  obtain  a  flexural  shear  wall  corresponding  to  the 
braced  truss;  instead  several  analyses  of  the  structure  were  performed  by 
varying  the  flexural  stiffness  of  the  shear  wall.  In  all  cases,  the 
ratio  of  the  wall  stiffness  to  the  column  stiffness  was  constant  in  each 
story  of  the  structure. 

The  plastic  moment  capacity  at  the  bottom  of  the  wall  was 
selected  so  to  bear  a  reasonable  relationship  to  the  wall  stiffness.  The 
plastic  moment  capacity  of  the  shear  wall  in  the  stories  above  the  base 
was  reduced  linearly  so  that  the  ratio  of  wall  strength  to  stiffness  would 
be  constant  for  each  story.  A  bilinear  moment-curvature  relationship  was 
used  for  the  wall  with  the  slope  of  the  strain-hardening  branch  being 
one-sixth  of  the  elastic  slope.  The  modulus  of  elasticity  of  the  frame 
was  30,000  kip/in^  and  for  the  shear  wall  3,000  kip/in^  .  A  constant 
8'-0"  wall  width  was  maintained  throughout  the  height  of  the  structure. 
High  values  of  the  rotational  spring  constants  were  used  at  the  base  of 
the  column  and  shear  wall  in  order  to  approach  a  fixed  end  condition. 
Concentrated  lateral  loads  were  assumed  to  act  at  each  floor  level.  The 
lateral  load  at  the  top  floor  was  one-half  that  at  the  other  levels.  The 
vertical  load  was  held  constant  during  the  analysis  and  the  lateral  loads 
were  increased  monotonically . 

First  and  second  order  (with  P-A  effect)  inelastic  analyses 
were  performed  for  the  24 -story  structure.  FIGURE  5.5  shows  graphs  of 
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the  lateral  force  at  the  top  of  the  frame,  H  ,  versus  the  sway  rotation 

of  the  top  story  (A  0/l  -  A  00)  /H  0/l  .  The  frame  has  a  ratio  of  wall 
r  J  w  24  w  23  s  24 

stiffness  to  column  stiffness  of  50  .  The  upper  curve  has  been  obtained 

from  a  first  order  inelastic  analysis.  The  first  hinge  in  the  frame  is 
detected  at  point  'a'  .  The  shear  wall  yields  first  at  the  base  as 

shown  by  point  ' b 1  on  the  graph.  The  structure  is  essentially,  a  'weak 
beam-strong  column'  type  and  at  point  ' c '  hinges  have  formed  at  the  ends 
of  all  the  beams.  The  only  column  hinge  detected  between  points  ' b ' 
and  ' c '  occurs  at  the  top  of  the  column  in  the  24th  story.  Since  a 
'strain-hardening'  type  of  moment-curvature  relationship  has  been  assumed 
for  the  wall,  the  wall  continues  to  accept  additional  load  .  To 
demonstrate  the  P-A  effect,  an  analysis  represented  by  the  lower  curve 
in  FIGURE  5.5  was  performed.  In  this  case,  the  structure  was  analyzed 
with  a  reduced  plastic  moment  capacity  (TABLE  5.4)  for  the  columns.  At 
point  'a'  on  the  lower  curve,  the  first  hinge  in  the  structure  was 
detected.  Up  to  point  ' b '  the  rate  of  decrease  of  the  stiffness  of 
the  structure  is  moderate,  however,  beyond  point  ' b '  the  deflection  of 
the  structure  increases  rapidly.  At  point  ' c '  most  of  the  wall  is 
inelastic.  For  the  next  increment  of  lateral  load,  the  deformations 
become  so  large  that  the  system  does  not  converge.  The  load  corresponding 
to  ' c '  has  been  taken  as  the  ultimate  load  carrying  capacity  of  the 
structure. 

In  FIGURE  5.5,  the  effect  of  the  P-A  moment  on  the  structure 
is  illustrated.  The  difference  in  load  carrying  capacity  predicted  by 
the  two  analyses  is  primarily  due  to  the  P-A  effect,  since  only  one 
hinge  forms  in  the  column  between  points  ' b '  and  ' c '  according  to  the 


47 


first  order  analysis  and  hinges  did  not  form  in  the  columns  according  to 
the  second  order  analysis. 

FIGURE  5.6  consists  of  plots  of  the  story  number,  N  ,  versus 
the  lateral  deflection,  A  ,  showing  the  deflected  shape  of  the  structure 
as  the  lateral  load  is  incremented.  The  curves  are  obtained  from  the 
second  order  analysis  in  FIGURE  5.5  The  values  of  the  applied  lateral 
load  at  the  top  of  the  structure  are  used  to  distinguish  the  curves.  The 
deflected  shape  for  loads  of  5.40  kips  and  5.60  kips  correspond  to 
steps  'b1  and  ' c '  respectively  in  FIGURE  5.5.  These  curves  emphasize 
the  rapid  increase  in  deflections  which  occurs  as  the  wall  enters  the 
inelastic  range. 

The  shear  resisted  by  the  frame  in  the  upper  stories  is  consid¬ 
erably  higher  than  the  applied  shear.  FIGURE  5.7  plots  the  top  level 
lateral  force,  H  ,  versus  the  percentage  of  the  story  shear  resisted  by 
the  frame,  V^/Vs  .  The  curve  is  obtained  from  the  second  order 
analysis.  In  FIGURE  5.7,  the  relationship  is  almost  straight  up  to  point 
'a'  since  few  hinges  have  formed  in  the  frame.  From  'a'  to  'c'  the 
percentage  shear  resisted  by  the  frame  increases  with  an  increase  in  the 
lateral  load.  Between  points  ' b 1  and  'c'  the  increase  is  quite  rapid. 
Points  1 b '  and  'c'  in  FIGURE  5.7  correspond  to  the  same  points,  in 
FIGURE  5.5. 

The  effect  of  varying  the  wall  stiffness  (keeping  the  wall 
width  constant)  was  also  studied.  FIGURE  5.8  is  a  plot  of  top  level 
lateral  force,  H  ,  versus  the  top  level  lateral  deflection,  A  .  In 
all  cases,  the  P-A  effect  (corresponding  to  the  total  vertical  loads 
listed  in  TABLE  5.4)  was  included  and  the  reduced  plastic  moment 
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capacities  were  used  for  the  columns. 

The  curves  in  FIGURE  5.8  represents  ratios  of  the  wall  to  column 
stiffness  ranging  from  5  to  50  .  This  ratio  was  constant  for  all 
stories.  The  load  carying  capacity  of  the  structure  decreases  as  the 
ratio  of  the  stiffness  of  the  shear  wall  to  the  column  is  decreased. 

The  structure  having  S^/ Sc  =  50  was  reanalyzed  with  zero  wall 
width.  The  result  is  shown  in  FIGURE  5.8  as  the  dashed  curve  .  The 
analysis  assuming  zero  wall  width  neglects  the  extra  restraining  moment 
on  the  wall.  This  moment  is  the  result  of  the  shear  at  the  wall  end  of 
the  right  beam  acting  through  half  the  wall  width.  For  the  example 
considered  the  analysis  predicts  a  load  carrying  capacity  below  that  of 
the  actual  structure. 

It  can  be  observed  from  FIGURE  5.8  that  the  difference  in 
behavior  due  to  the  variation  in  wall  to  column  stiffness  is  considerable. 
In  all  cases,  the  analysis  of  the  shear  wall-frame  structure  did  not 
converge  beyond  point  ' c '  .  In  the  analysis,  the  loads  corresponding 

to  points  1 c 1  have  been  taken  as  the  ultimate  load  carrying  capacities 
of  the  structures.  Due  to  the  procedure  used  in  the  analysis,  the 
unloading  branch  of  the  load-deflection  curve  for  the  structure  can  not 


be  obtained. 


49 


STORY  NO.  VERTICAL  MOMENT  OF  INERTIA  (IN4) 

OR  LOADS 

FLOOR  NO. 

(kips)  BEAM  COLUMN  WALL 


1 

6. 

25 

X 

104 

4.66 

X 

5 

10° 

9.50 

X 

5 

10° 

5.02 

X 

10 

2 

5. 

75 
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104 

4.66 

X 

5 

103 

5.96 

X 

105 

5.02 

X 

10 

3 

5. 

25 
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4.66 

X 

5 

10° 

5.70 

X 

105 

5.02 
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10 

4 
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4.58 
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5.38 
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5.02 

X 

10 
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X 
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4.50 

X 

5 

10° 

4.45 

X 
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5.02 

X 

10 

8 

2. 

93 

X 

104 

4.50 

X 

5 

10° 

4.16 

X 

105 

5.02 

X 

10 

9 

2. 

48 

X 

104 

4.50 

X 

105 

3.84 

X 
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5.02 

X 

10 

10 

2. 

05 

X 
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4.42 

X 
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3.54 

X 
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5.02 

X 

10 

11 
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3.24 
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5.02 
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10 

12 
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4.42 

X 
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5.02 
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13 

8. 

00 
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X 
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10° 

5.02 

X 

10 

14 
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X 
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X 
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5.02 

X 

10 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 
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TABLE  5.1 


LOADS  AND  MEMBER  PROPERTIES  (FOURTEEN-STORY  STRUCTURE) 


* 

MOMENT  OF  INERTIA  OF  WALL  (STORY)  PERCENTAGE  STORY  SHEAR 
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LEVEL 

MOMENT  OF 

INERTIA  (IN4) 

PLASTIC  MOMENT 

CAPACITY  (KIP-IN) 

BEAM  AB 

BEAM  BC 

BEAM  AB 

BEAM  BC 

1 

888 

1170 

5404 

4920 

2 

888 

1170 

5404 

4920 

3 

888 

1170 

5404 

4920 

4 

888 

1170 

5404 

4920 

5 

888 

984 

5404 

4560 

6 

888 

984 

5404 

4560 

7 

888 

984 

5404 

4560 

8 

888 

984 

5404 

4560 

9 

888 

890 

5404 

4020 

10 

888 

890 

5404 

4020 

11 

888 

890 

5404 

4020 

12 

888 

890 . 

5404 

4020 

13 

888 

890 

5404 

4020 

14 

888 

890 

5404 

4020 

15 

888 

890 

5404 

4020 

16 

888 

890 

5404 

4020 

17 

888 

890 

5404. 

4020 

18 

888 

890 

5404 

4020 

19 

888 

890 

5404 

4020 

20 

888 

890 . 

5404 

4020 

21 

888 

890 

5404 

4020 

22 

888 

890 

5404 

4020 

23 

888 

890 

5404 

4020 

24 

630 

584 

3900 

2960 

TABLE  5.3  BEAM  PROPERTIES  (TWENTY  FOUR-STORY  STRUCTURE) 


1 

OR 

STO 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 
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COLUMN  PROPERTIES 


MOMENT  OF  PLASTIC  VERTICAL 

INERTIA  MOMENT 

A  CAPACITY  LOAD 

( IN  J  (K-IN)  (KIPS) 


REDUCED  PLASTIC 
MOMENT-CAPACITY 
OF  COLUMN 

(KIP-IN) 


14814 

73000 

12680 

27300 

14814 

73000 

12330 

27300 

13278 

66100 

11600 

25400 

13278 

66100 

11000 

25400 

11452 

58000 

10540 

21900 

11452 

58000 

10000 

21900 

9716 

50600 

9440 

18400 

9716 

50600 

8930 

18400 

8306 

43600 

8400 

15750 

8306 

43600 

7860 

15750 

7098 

37800 

7320 

13500 

7098 

37800 

6790 

13500 

5850 

31600 

6260 

11500 

5850 

31600 

5720 

11500 

4850 

26600 

5180 

10800 

4850 

26600 

4650 

10800 

4130 

22900 

4120 

10300 

4130 

22900 

3590 

10300 

2970 

16700 

3050 

6960 

2970 

16700 

2520 

6960 

1620 

10500 

1980 

5580 

1620 

10500 

1450 

5580 

1100 

6060 

1050 

5260 

1100 

6060 

370 

5260 

TABLE  5.4  VERTICAL  LOADS  AND  COLUMN  PROPERTIES 
(TWENTY  FOUR-STORY  STURCTURE) 
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FIGURE  5.1  FOURTEEN-STURY  STRUCTURE 
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FIGURE  5.3  DEFLECTED  SHAPE  OF  FOURTEEN-STORY  STRUCTURE 
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FIGURE  5.4  TWENTY-FOUR-STORY  STRUCTURE 
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FIGURE  5.5  TOP  LEVEL  LOAD  -  TOP  STORY  SWAY  CURVE  (WALL  TO  COLUMN  STIFFNESS  -  50) 


FIGURE  5.6  DEFLECTED  SHAPE  OF  TWENTY-FOUR-STORY  STRUCTURE 


(WALL  TO  COLUMN  STIFFNESS  -  50) 


FIGURE  5.7  VARIATION  OF  TOP  STORY  FRAME  SHEAR 


A  ( INCHES) 

FIGURE  5.8  TOP  LEVEL  LOAD-DEFLECTION  CURVE  (TWENTY-FOUR-STORY  STRUCTURE) 


CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

A  method  of  second  order  inelastic  analysis  of  combined  shear 
wall-frame  structures  has  been  presented  herein.  A  computer  program  was 
developed  to  facilitate  the  computations  involved. 

In  this  analysis  the  actual  structure  has  been  replaced  by  a 
simplified  structural  model.  The  problems  involved  in  determining  a 
practicable  model  are  discussed  with  reference  to  previous  work.  A 
conservative  method  has  been  used  for  lumping  the  actual  structure  into 
the  model.  However,  the  justification  for  using  this  method  for  various 
types  of  shear  wall-rigid  frame  structures  is  not  yet  fully  developed. 

The  analysis  of  the  model  is  discussed  in  detail.  The  model 
was  analyzed  taking  into  account  the  formation  of  plastic  hinges  in  the 
frame,  the  inelastic  action  of  the  shear  wall  and  the  secondary  moments 
caused  by  the  P-A  effect.  The  method  of  analysis  consists  basically 
of  three  iterative  processes:  the  first  determines  the  joint  rotations 
of  the  displaced  frame  system;  the  second  determines  the  distribution  of 
lateral  shear  between  the  shear  wall  and  the  frame;  and  the  third  takes 
into  account  the  P-A  effect. 

In  these  processes,  the  inelastic  action  of  the  wall  is  included 
by  modifying  the  stiffness  of  the  yielded  segments  of  the  wall  so  as  to  be 
consistent  with  the  assumed  moment-curvature  relationship.  The  formation 
of  plastic  hinges  in  the  frame  is  accounted  for  by  modifying  the  joint 
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rotation  equations  at  the  pertinent  locations.  After  the  formation  of  a 
plastic  hinge  in  the  frame,  the  analysis  is  continued  by  assuming  that  an 
imaginary  hinge  is  inserted  at  the  cross-section;  this  hinge  is  subjected 
to  opposing  moments  equal  to  the  plastic  moment  capacity  of  the  section. 
The  P-A  effect  has  been  simulated  by  subjecting  the  analytical  model  to 
an  equivalent  lateral  force  system,  which  is  increased  to  account  for  the 
presence  of  the  secondary  moments. 

A  printout  of  the  computer  program  is  given  in  APPENDIX  C.  The 
results  of  the  analyses  of  two  example  frames  have  been  presented  and 
discussed. 

It  is  concluded  that  the  analytical  model  yields  comparable 

results  (in  the  elastic  range)  to  those  obtained  by  using  the  method 

(2) 

developed  previously  v  .  The  model  is  relatively  simple  and  gives  a  good 
indication  of  the  behavior  of  the  structure.  The  rationality  of  the 
lumping  procedure  used  and  the  behavior  of  the  structure  in  the  inelastic 
range  are  still  to  be  checked  by  more  rigorous  analysis  and  tests. 

The  analyses  performed  have  clearly  shown  the  necessity  of 
considering  the  interaction  between  the  frame  and  the  shear  wall.  This 
interaction  must  be  evaluated  if  the  actual  force  distribution  in  the 
structure  is  to  be  estimated.  Generally  in  tall  structures,  the  shears 
developed  in  the  top  stories  are  much  higher  than  the  applied  story 
shears.  Neglect  of  this  factor  could  lead  to  excessive  deflections  and 
possible  failure  of  the  affected  frame  members. 

The  P-A  effect  caused  aconsiderable  reduction  in  load 
carrying  capacity  for  the  more  flexible  structure  analyzed.  The  secondary 
moments  produced  by  the  P-A  effect,  combined  with  the  inelastic  action 
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of  the  wall  and  frame  resulted  in  reduced  load-carrying  capacities.  The 
frames  analyzed  did  not  form  a  mechanism  at  the  ultimate  load  but  instead 
failure  was  due  to  overall  instability  of  the  structure. 

The  vertical  deflections  of  the  beam-to-col umn  joints  due  to 
axial  column  shortening  have  not  been  included  in  this  analysis.  This 
effect  may  influence  the  deformations  of  the  structure  resulting  in  an 
additional  reduction  in  the  ultimate  load  carrying  capacity. 

Although  the  analysis  has  been  checked  by  manual  computations 
on  simple  structures,  it  is  recommended  that  the  rationality  of  the 
assumptions  made  be  verified  by  the  following  three  steps:  first,  the 
subassemblage  analysis  should  be  checked  against  a  rigorous  analysis 
performed  by  using  an  independent  program,  secondly,  the  actual  structure 
should  be  analyzed  rigorously  and  the  results  compared  with  those  of  a 
subassemblage  analysis,  thirdly,  the  predicted  behavior  should  be  checked 
by  performing  full  scale  tests. 

Due  to  the  method  of  analysis  used,  it  is  not  possible  to  obtain 
the  unloading  branch  of  the  load-deflection  curve.  In  order  to  understand 
thoroughly  the  behavior  of  the  structure,  the  unloading  branch  must  be 


determined. 
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D. .  . 
W  l 

E 

F  , 


FW  i 


H. 

Hi  i 


H 


S  1 


K 


12 


Kr 


K 


BE 


K, 


K 


CE 


12 


M 


M 


pc 

BF  i 


NOMENCLATURE 

Width  of  shear  wall  system  at  the  i  th  floor 
Youngs  modulus  of  elasticity 

Applied  force  on  the  wall  at  the  i  th  floor  from  the 
frame  analysis 

Portion  of  the  applied  lateral  load  resisted  by  the 
wall  at  the  i  th  floor 

Applied  lateral  load  on  the  structure  at  the  i  th  floor 
Additional  lateral  load  to  be  applied  at  the  i  th  floor 
in  the  first  cycle  to  simulate  the  P-A  effect  on  the 
structure 

Story  height  in  the  i  th  story 

End  moment  of  the  beam  (12)  at  the  end  (1)  for  a  unit 
rotation  at  the  end  (2) 

Stiffness  of  one  beam 
Equivalent  beam  restraint 
Stiffness  of  one  column 
Equivalent  column  stiffness 
Length  of  beam  (12) 

Plastic  moment  capacity  of  shear  wall 

Plastic  moment  capacity  of  column  reduced  for  axial  load 

Moment  at  the  column  end  of  the  left  hand  beam  at  the 


i  th  floor 
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MgW  .  Moment  at  the  column  end  of  the  right  hand  beam  at  the 

i  th  floor 

Mc  .  Moment  at  the  top  of  the  column  of  the  i  th  story 

Mc  .+i  Moment  at  the  bottom  of  the  column  of  the  i+1  th  story 

M. ,  •  Total  moment  on  the  wall  at  the  i  th  floor  from  the 

frame  analysis 

M^b  .  Moment  at  the  wall  end  of  the  right  hand  beam  at  the 

i  th  floor 

MPgp  .  Plastic  moment  capacity  of  the  left  hand  beam  at  the 

i  th  floor 

MPg^j  .  Plastic  moment  capacity  of  the  right  hand  beam  at  the 

i  th  floor 

MP  .  Plastic  moment  capacity  (reduced  for  axial  load)  of  the 

v*  I 

column  of  the  i  th  story 

P.  Total  vertical  load  on  the  structure  above  the  i  th  floor 

R  .  Chord  rotation  of  the  column  of  the  i  th  story 

R'  .  Vertical  chord  rotation  of  the  right  hand  beam  connecting 

the  column  to  the  wall  at  the  i  th  floor 

SD[-  •  Stiffness  (EI/L)  of  the  left  hand  beam  at  the  i  th  floor 

Dr  1 

S  .  Stiffness  (EI/L)  of  the  column  of  the  i  th  story 

c  i 

V'  .  Extra  shear  in  the  structure  at  the  i  th  story  due  to 

P-A  effect 

A-j  .  Final  deflection  of  the  structure  under  the  lateral  load 

system  H.  at  the  i  th  floor.  The  subscript  1  refers 
to  the  first  cycle  in  the  P-A  iteration  process 


- 


. 
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A2  i 

Final  deflection  of  the  structure  under  the  lateral  load 

A  (n) 

AWI 

system  H.  +  H-j 1  .  at  the  i  th  floor 

Lateral  displacement  of  the  wall  at  the  i  th  floor  in  the 

A'  (n) 

AWI 

n  th  cycle  of  the  wall -frame  iteration  process 

Lateral  displacement  of  the  wall  at  the  i  th  floor  after 

applying  forced  convergence  formula  in  the  n  th  cycle 

of  the  wall-frame  iteration  process 

0  . 

1 

Rotation  of  the  frame  joint  at  the  i  th  floor 

9F  i 

Rotation  at  the  left  end  of  the  left  hand  beam  at  the 

i  th  floor 

0 

wW  i 

Slope  of  the  wall  at  the  i  th  floor  in  the  n  th  cycle  of 

0' 

H  W  i 

the  wall -frame  iteration  process 

Slope  of  the  wall  at  the  i  th  floor  after  applying  forced 

convergence  formula  in  the  n  th  cycle  of  the  wall -frame 

iteration  process 

9WB  i 

Rotation  at  the  wall  end  of  the  right  hand  beam  at  the 

i  th  floor 

v  pc 

Curvature  corresponding  to  M  assuming  ideally  elastic 

pc 

,  (n) 

5WI 

behavior  (0  =  M  /El ) 

p  p 

Vertical  deflection  of  the  beam-to-column  connection  at 

the  i  th  floor,  in  the  n  th  cycle  of  the  wall -frame 

A' 

6  WI 

iteration  process 

Vertical  deflection  of  the  beam-to-wall  connection  at  the 

i  th  floor  in  the  n  th  cycle  of  the  wall -frame  iteration 

process 


' 
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BAKA 

BMOM 

BMOMF 

CON 

CDEFF 

CROTF 

DEF 

DEFF 

DMBF 

DMBW 

DMCB 

DMCT 

DMFB 

DMWB 


NOMENCLATURE  FOR  FORTRAN  IV  PROGRAM 

A  SUBROUTINE  SUBPROGRAM 

BOTTOM  MOMENT  IN  A  SEGMENT  OF  THE  SHEAR  WALL  (KIP-IN) 
BOTTOM  MOMENT  IN  THE  BOTTOM  OF  THE  LOWEST  SEGMENT  IN 
A  STORY  (KIP-IN) 

CONVERGENCE  LIMIT  FOR  THE  DEFLECTION  AND  ROTATION 
FLOOR  LEVEL  DEFLECTION  OF  SHEAR  WALL  AFTER  APPLYING 
CONVERGENCE  FORMULA 

FLOOR  LEVEL  ROTATION  OF  SHEAR  WALL  AFTER  APPLYING 

CONVERGENCE  FORMULA 

DEFLECTION  OF  WALL  AT  A  SECTION  (IN) 

DEFLECTION  OF  THE  WALL  AT  THE  FLOOR  LEVEL  (IN) 

MOMENT  AT  THE  COLUMN  END  OF  THE  LEFT  HAND  BEAM 
COMPUTED  IN  THE  PREVIOUS  LOADING  CONDITION  (KIP-IN) 
MOMENT  AT  THE  COLUMN  END  OF  THE  RIGHT  HAND  BEAM  COMPUTED 
IN  THE  PREVIOUS  LOADING  CONDITION  (KIP-IN) 

MOMENT  AT  THE  BOTTOM  END  OF  COLUMN  IN  A  PARTICULAR  STORY 
COMPUTED  IN  THE  PREVIOUS  LOADING  CONDITION  (KIP-IN) 
MOMENT  AT  THE  TOP  END  OF  A  COLUMN  IN  A  PARTICULAR  STORY 
COMPUTED  IN  THE  PREVIOUS  LOADING  CONDITION  (KIP-IN) 
MOMENT  AT  THE  LEFT  END  OF  THE  LEFT  HAND  BEAM  COMPUTED 
IN  THE  PREVIOUS  LOADING  CONDITION  (KIP-IN) 

MOMENT  AT  THE  WALL  END  OF  THE  RIGHT  HAND  BEAM  COMPUTED 
IN  THE  PREVIOUS  LOADING  CONDITION  (KIP-IN) 


•  Md 
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DW 

EF 

EW 

F 

FD 

FRAME 

FOC 

FW 

HBF 

HBW 

HLI 

HLIR 

HS 

HSF 

IK  ,  IX 
IVL 

KB 

KC 

L 

MAX 

MI 

MIBF 


WIDTH  OF  SHEAR  WALL  (IN) 

YOUNGS  MODULUS  OF  ELASTICITY  OF  THE  FRAME  (KIP/IN2) 

YOUNGS  MODULUS  OF  ELASTICITY  OF  THE  WALL  (KIP/IN2) 

APPLIED  FORCE  ON  THE  WALL  FROM  FRAME  ANALYSIS  (KIP) 
HORIZONTAL  DESIGN  LOAD  ACTING  AT  FLOOR  LEVEL  (KIP) 

A  SUBROUTINE  SUBPROGRAM 

FRAME  FORCE  AT  THE  FLOOR  LEVEL  (KIP) 

FINAL  FORCE  ON  WALL  AT  EACH  FLOOR  LEVEL  (KIP) 

LENGTH  OF  LEFT  HAND  BEAM  (IN) 

LENGTH  OF  RIGHT  HAND  BEAM  (IN) 

PERCENTAGE  OF  ORIGINAL  HORIZONTAL  LOAD  TO  BE  INCREASED 
IN  THE  ELASTIC  RANGE 

PERCENTAGE  OF  ORIGINAL  HORIZONTAL  LOAD  TO  BE  INCREASED 
IN  THE  INELASTIC  RANGE 
STORY  HEIGHT  (IN) 

HEIGHT  OF  FLOOR  LEVEL  FROM  BASE  (IN) 

DUMMY  CONSTANTS 

DUMMY  CONSTANTS  FOR  STOPPING  THE  PROGRAM  IF  THE  DEFORMATION 
EXCEEDS  CERTAIN  SPECIFIED  LIMIT 

SPRING  CONSTANT  AT  THE  BASE  OF  THE  SHEAR  WALL  (KIP-IN/RAD) 
SPRING  CONSTANT  AT  THE  BASE  OF  THE  SHEAR  WALL  (KIP-IN/RAD) 
PRODUCT  OF  NUMBER  OF  STORY  AND  NUMBER  OF  SEGMENTS  IN  A  STORY 
MAXIMUM  NO.  OF  CYCLE  TO  BE  PERFORMED  FOR  ANY  ITERATION 
PROCESS 

MOMENT  OF  INERTIA  OF  SHEAR  WALL  (IN4) 

MOMENT  OF  INERTIA  OF  LEFT  HAND  BEAM  (IN4) 
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MIBW 

MIC 

Mil 

MM 

MMCB 

MMCT 

MMMM 

MOMB 

MOMBF 

MOMBW 

MOMFB 

MOMP 

MOMW 

MOMWB 

MPC 

MPF 

MPSW 

MPW 

ND 

NI 


MOMENT  OF  INERTIA  OF  RIGHT  HAND  BEAM  (IN4) 

MOMENT  OF  INERTIA  OF  COLUMN  (IN4) 

MOMENT  OF  THE  SEGMENT  OF  THE  SHEAR  WALL  (IN4) 

NO.  OF  PROBLEMS  TO  BE  SOLVED 

MOMENT  AT  THE  BOTTOM  OF  A  COLUMN  IN  A  STORY  (KIP-IN) 

MOMENT  AT  THE  TOP  OF  A  COLUMN  IN  A  STORY  (KIP-IN) 

COUNTER  FOR  ITERATION  WITH  THE  EFFECT  OF  AXIAL  LOAD 
MOMENT  AT  THE  BASE  OF  THE  SHEAR  WALL  (KIP-IN) 

MOMENT  AT  THE  COLUMN  END  OF  THE  LEFT  HAND  BEAM  IN  THE 

CYCLE  UNDER  CONSIDERATION  (KIP-IN) 

MOMENT  AT  THE  COLUMN  END  OF  THE  RIGHT  HAND  BEAM  IN  THE 
CYCLE  UNDER  CONSIDERATION  (KIP-IN) 

MOMENT  AT  THE  LEFT  END  OF  THE  LEFT  HAND  BEAM  IN  THE 
CYCLE  UNDER  CONSIDERATION  (KIP-IN) 

MOMENT  IN  A  STORY  DUE  TO  AXIAL  LOAD  (KIP-IN) 

TOTAL  MOMENT  AT  FLOOR  LEVEL  ON  THE  SHEAR  WALL  DUE  TO  END 
MOMENT  AND  SHEAR  FROM  THE  RIGHT  HAND  BEAM  (KIP-IN) 

MOMENT  AT  THE  RIGHT  END  OF  THE  RIGHT  HAND  BEAM  IN  THE 
CYCLE  UNDER  CONSIDERATION  (KIP-IN) 

PLASTIC  MOMENT  CAPACITY  OF  COLUMN  IN  A  PARTICULAR 
STORY  (KIP-IN) 

PLASTIC  MOMENT  CAPACITY  OF  THE  LEFT  HAND  BEAM  (KIP-IN) 
PLASTIC  MOMENT  CAPACITY  OF  THE  SHEAR  WALL  (KIP-IN) 
PLASTIC  MOMENT  CAPACITY  OF  THE  RIGHT  HAND  BEAM  (KIP-IN) 
NUMBER  OF  DIVISION  IN  A  STORY 
NUMBER  OF  TIMES  HORIZONTAL  LOAD  TO  BE  INCREMENTED 


(■Ml -9$)  M0T2  A  MI  IlfWJOO  A  ft  ^0-  3HT  TA  0I3WK 
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NM 

LENGTH  OF  A  SEGMENT  OF  A  WALL  IN  A  STORY  (IN) 

NNMM 

SEGMENT  HEIGHT  FROM  BASE  OF  WALL  (IN) 

NS 

NUMBER  OF  STORIES 

ROB 

SWAY  ROTATION  OF  RIGHT  HAND  BEAM 

ROFA 

A  SUBROUTINE  SUBPROGRAM 

ROS 

STORY  ROTATION  OF  FRAME 

RMWP 

RATIO  OF  MOMENT  TO  PLASTIC  MOMENT  CAPACITY  OF 

WALL 

RPWP 

RATIO  OF  CURVATURE  TO  CURVATURE  CORRESPONDING 

TO  RMWP 

ROT 

ROTATION  OF  WALL  IN  A  SEGMENT  (RAD) 

ROTF 

ROTATION  OF  WALL  AT  THE  FLOOR  LEVEL  (RAD) 

ROTFF 

JOINT  ROTATION  OF  FRAME  (RAD) 

ROTO 

JOINT  ROTATION  AT  THE  BASE  OF  COLUMN  (RAD) 

SBF 

STIFFNESS  OF  LEFT  HAND  BEAM  (KIP-IN) 

SBW 

STIFFNESS  OF  RIGHT  HAND  BEAM  (KIP-IN) 

SC 

STIFFNESS  OF  COLUMN  (KIP-IN) 

SHEARR 

SHEAR  AT  THE  ENDS  OF  LEFT  HAND  BEAM  (KIP) 

SHEARW 

SHEAR  AT  THE  ENDS  OF  RIGHT  HAND  BEAM  (KIP) 

SHEC 

SHEAR  IN  COLUMN  IN  A  STORY  (KIP) 

SR 

A  SUBROUTINE  SUBPROGRAM 

SR  1 

A  SUBROUTINE  SUBPROGRAM 

SR  3 

A  SUBROUTINE  SUBPROGRAM 

SSB 

SLOPE  OF  THE  SECOND  BRANCH  OF  SHEAR  WALL  MOMENT-CURVATURE 

DIAGRAM 

TMOMF 

TOP  MOMENT  IN  THE  TOP  OF  THE  TOPMOST  SEGMENT 

(KIP-IN) 

IN  A  STORY 
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TMOM 

VDEF 

PPS 


TOP  MOMENT  IN  A  SEGMENT  OF  THE  SHEAR  WALL  (KIP-IN) 
VERTICAL  UPLIFT  OF  THE  BEAM  CONNECTED  WITH  THE  WALL  (IN) 
TOTAL  AXIAL  LOAD  IN  STRUCTURE  AT  A  FLOOR  LEVEL  (KIP) 
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APPENDIX  A 

DERIVATION  OF  JOINT  ROTATION  EQUATION 

The  frame  system  has  been  forced  into  the  deformed  shape  of 
the  wall  shown  in  FIGURE  4.5.  Applying  the  slope-deflection  equations, 
the  moment  MgW  .  ,  at  the  end  B.  of  the  right  beam  is  given  by 


M 


BW  i 


=  S 


BW  i 


6i  +  2  0UB  i 


6  R  ‘  i ) 


where  .  =  EI/L  is  the  stiffness  of  the  right  beam  (B..  W. )  and 

R1.  is  the  vertical  chord  rotation  of  the  right  beam  connecting  the 
column  to  the  wall.  Also  in  the  above  equation  0.  is  the  rotation 
of  the  frame  joint  and  Q,m  .  is  the  rotation  at  the  wall  end  of  the 

Wd  1 

right  beam. 

Considering  the  left  beam  B..F..  and  applying  the  slope 

deflection  equation,  assuming  that  the  joint  rotation  at  the  two  ends 

B.  and  F.  are  equal,  the  moment,  MDC  .  ,  at  the  end  F.  is 

i  i  M  BF  l  i 

given  by 


M 


BF  i 


=  6.  S 


BF  i 


where  S 


BF  i 


EI/L  is  the  stiffness  of  the  left  hand  beam, 


B.  F. 

l  l 
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The  column  B.  B^-j  between  the  i  th  and  the  i+l  th  story 

.  ,  ....  r  c  and  the  chord  rotation  of  this  column  is 

has  a  stiffness  of  S 

c  i  +1 

R.^n  .  The  moment,  M  . , n  ,  at  the  bottom  end  B.  of  the  column 
i+l  c  i+ I  1 

Bi  B.+1  is  given  by 

M„  =  S  • , -i  (4  0.  +  2  0.,,  -  6  R .  ) 

c  i+l  c  i+l  v  1  i+l  i+l 

where  0.,,  is  the  rotation  of  the  frame  joint  at  the  i+l  th  level, 
i  +  l 

The  moment,  M  .  at  the  top  end,  B.  ,  of  the  column 

Cl  5  r  i 

B.  B._.|  is  given  by: 


M 


c  1 


=  S 


c  1 


(4  0. 


+  2  9.  , 
i-l 


6  R.) 


where  S  .  =  EI/L  is  the  stiffness  of  the  column  B.  B.  and  R. 

Cl  1  1-1  1 

is  the  chord  rotation  of  the  column.  0.  i  is  the  rotation  of  the 
frame  joint  at  the  i-l  th  level. 

Equating  the  sum  of  the  joint  moments  to  zero 


M 


BW  i 


+  M 


BF  i 


+  M. 


c  1 


+  M 


c  i+l 


=  0 


Substituting  the  values  of  the  end  moments  from  the  above 


equations  the  joint  equilibrium  equation  becomes 


* 
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SBW  i  (4  9i  +  2  0WB  1  ‘  6  RV  +  6  SBF  i  ei  + 

sci  <40i+20i-l  -  6  Ri)  +  Sc  1+1 
(4  0.  +  2  6.+1  -  6  Ri+1)  -  0 

and  after  simplification: 


Si  ■  {SBW  i  <6R'  i  -  20WB  1>  +Sc  i  <6  Ri  -  20i-l> 


+  Sc  1+1  (6  Ri+1  -  20i+l>}  /  (4  SBW  i  +  4  Sc  i  + 


4  Sc  i+1  +  6  SBF  i} 


(A-l) 


form 


Equation  (A-l)  can  be  written  in  the  following  simplified 


q  _A  +  B  +  C  +  D 

”  A  I 


’i  AT  +  B1  +  C'  +  D1 


(A-2) 


Where  in  the  numerator, 


A  ~  S BW  i  Rl  i  ■  2  0WB  i ^ 

B  =  sc  t  (6  R.  -  2  9  i_1) 

C  =  Sc  i+1  Ri+1  "  2  0  i+1 ^ 


■ 
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and 


D  =  0 


And  in  the  denominator, 


A1  =  4  S 


BW  i 


C'  =  4  S 


c  i+1 


and 


D'  =  6  S 


BF  i 


Equations  (A-l)  and  (A-2)  are  valid  for  the  elastic  analysis 


only.  When  a  plastic  hinge  forms  at  any  one  of  the  potential  hinge 
locations  (FIGURE  4.7),  the  joint  rotation  equation  must  be  modified. 
Table  A.l  lists  the  substitutions  to  be  made  in  Equation  (A-2)  for  the 
formation  of  hinges  in  the  various  members,  so  that  the  joint 
rotation  equation  will  conform  to  the  particular  hinge  pattern 
considered. 


In  Table  A.l,  MPgp  .  ,  MP^  .  are  the  plastic  moment 

capacities  of  the  beams  B.  F.  and  B.  W.  ,  and  mp  .is  the 
r  11  11  c  i 


plastic  moment  capacity  (reduced  for  axial  load)  of  the  column. 
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TABLE  A.  1 

MODIFICATION  OF  ELASTIC  SLOPE-DEFLECTION 
EQUATIONS  FOR  JOINT  EQUILIBRIUM 


Member 


B.  W. 
i  i 


Hinge  Location  Substitution  to  be  Made 

In  Equation  (A-2) 


W. 


A  "  3  SBW  i  R' i  "  MPbw  i  1  2 


A'  =  3  S 


BW  i 


OR 


B.  ,  W. 

l  i 


A  =  MP 


BW  i 


A1  =  0 


B  =  3  Sc  j  •  Ri  ’  MPC  .  /  2 


B' 


3  S 


c  i 


B.  B 

i 


i-1 


B.  B  =  MP  . 

i  c  i 


OR 


B 


*  Bi-1 


B '  =  0 
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Member 

Hinge  Location  Substitution  to  be  Made 

In  Equation  (A-2) 

Bi  Bi+1 

Bi+1  c  =  3  sc  1+1  '  Ri+1  '  MPC  i+1  /  2 

C'  ■  3  Sc  i+1 

Bi  C  ■  MPc  i+1 

OR 

Bi+1  C'  =  0 

Bi  Fi 

Fi  D  ■  MPBF  i  '  2 

D1  3  SBF  . 

Bi  D  =  MPBF  i 

OR 

B.  ,  F.  D '  =  0 

i  i 

D 1  =  0 
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APPENDIX  B 

INCLUSION  OF  THE  P-A  EFFECT 

The  P-A  effect  plays  a  dominant  role  in  the  behavior  of 
flexible  structures.  The  vertical  load  on  the  structure,  P  ,  acting 
through  a  sidesway  displacement,  A  ,  produces  an  additional  overturn¬ 
ing  moment  commonly  known  as  the  P-A  moment.  Consider  a  single  story 
structure,  which  has  been  displaced  by  an  amount  A  and  is  subjected 
to  vertical  loads,  P  ,  as  shown  in  FIGURE  B.l.  To  balance  the  P-A 
moment  an  additional  shear  of  Pa/Hs  is  required  in  the  columns  of  the 
frame. 

In  this  analysis,  the  P-A  effect  has  been  taken  into 
account  by  analyzing  the  structure  under  an  equivalent  lateral  force. 

In  FIGURE  B.2a,  H.  is  the  applied  lateral  load  at  the  i  th  floor. 

Due  to  this  applied  lateral  force  system,  a  lateral  deflection  a-j  . 
will  be  produced. 

The  story  moment  between  the  i  th  and  i-1  th  floor  due  to 
the  P-A  effect  is  obtained  by  taking  moments  about  the  i  th  floor 
of  all  the  axial  loads  above  this  level.  Similarly,  by  taking 
moments  at  each  level,  the  other  story  moments  are  obtained.  The 
extra  shear  due  to  the  P-A  effect,  VI  ,  is  given  by: 

v  =  Mii  1  i-l) 

i  H  . 


■ 

■ 
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where  H$  .  is  the  story  height  and  P  is  the  total  vertical  load 
on  the  structure  above  the  i  th  floor. 

The  additional  lateral  load,  .  ,  simulating  the  P-A 

effect  can  then  be  determined  and  is  given  by 

H'  i  =  vi  -  v-!+i  =  Pi  ^  1  "  A]  i_1^  /  Hs  i 


Pi+1  ^A1  i+1  "  A1  i^ 


H 


s  i+1 


(B-l ) 


The  structure  is -now  reanalyzed  under  the  lateral  force 
system,  H.  +  .  (without  considering  vertical  loads).  In 

FIGURE  B.2b,  .  is  the  lateral  deflection  corresponding  to  the 

new  lateral  force  system.  The  deflections  .  etc.  are  then 

compared  with  the  deflections,  A-j  .  etc.  If  these  deflections  do  not 
agree  (within  a  specified  convergence  limit)  the  additional  horizontal 
force,  H£  •  ,  due  to  the  P-A  effect  for  the  next  cycle  is  given  by 

Hj  i  =  Pi^A2  i  '  A2  i-’P  .  Pi+1  W  i+1  "  A2  it.(B-2) 
Hs  i  Hs  i+l 

The  structure  is  again  analyzed  under  the  lateral  force 
system  H.  +  ^  .  The  process  is  continued  until  the  changes  in 

lateral  deflections  are  within  the  convergence  limit. 

The  total  vertical  load  at  each  floor  level  has  been  read 
into  the  program.  In  the  present  method  of  simulating  the  P-A  effect 


' 
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a  constant  vertical  load  has  been  used  throughout  the  analysis.  In  the 
actual  structure  the  girder  shears  will  introduce  tension  in  one  column 
and  compression  in  the  other.  This  will  not  change  the  total  vertical 
load  or  the  gross  P-A  effect. 
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X 


FIGURE  B . 1  P-A  EFFECT  -  SINGLE  STORY  FRAME  FIGURE  B.2  P-A  EFFECT  -  MULTI-STORY  STRUCTURE 


ooo  non  o  noon 
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FORTRAN  PROGRAM  FOR  THE  ANALYSIS  OF  SHEAR  WALL-FRAME  STRUCTURE 


DIMENSION  F( 99,30 ) ,HSF I  30 ) ,DW( 30} ,F0C ( 3 0 ) , SH E AR <  30 ) , ROT F ( 99 , 30 )  , 
1R0TFFI 99,30) , ROT 0 ( 99 ) , SHEC ( 30 ) , FD ( 3 0 ) , FW ( 30 ) , SHE ARW ( 3 0 ) , P S ( 30 ) , 
2SHEARR ( 30) ,CPHI ( 300) ,0(300) ,DIl 300) ,H{99, 30) , DELTA ( 99 , 30 ) , 

3FA( 30) , PPS( 30) ,CROTF (99,30) , CDEFF ( 99 , 30 ) ,  PPD<  30) , ZM { 300) , 

A  DEFF ( 99 , 30  ) ,  H8F { 30 ) , HBW ( 30 ) , HS ( 30 ) , SBF ( 30 ) , SC ( 3 0 )  , SB W ( 30 )  , 

5FI1( 30) , TMOM (300 ) ,BMOM( 300) , TMOMF (30) , BMOMF (30) ,DMFB( 30 ) , 

6DMBF(30) , DMBW ( 30  ),DMW8( 30) ,CSHEC(30) ,  DMCB(30) ,DMCT( 30) 

REAL  NNM( 300 ) ,NNMMC300) ,  M I  A ( 300 ) , M I  I { 300) , MMC T ( 30) ,MMCB( 30) , 

1KB, KC,  M0MW(30) ,  MOMWB ( 30 ) , MOMBW ( 30 ) , MOM BF ( 3 0 ) , MOMFB ( 30 ) 
2,MPF(30),MPw(30)  ,  MPC ( 30 ) , MRS ( 300 ) 

20  FORMAT ( 1H1) 

21  FORMAT ( 1HK ) 

22  FORMAT (1H  ) 

READ( 3, 1301  MM 

130  FORMAT! IX, 13 ) 

DO  400  J  J= 1 »  MM 

INPUT  STATEMENTS  FOR  DATA  BY  SUBROUTINE  SR. 

CALL  SR  ( L , NS , NO ,  EW  ,  KC  ,  KB  ,  FO  ,  N  I  ,  HS  ,  Dw  ,  P  S  ,  MPF , NNM , 

1  MPW ,MPC ,  SC, SBW,SBF, IVL,PPD,  HBF , HBW , 

2  RMWP,RPWP,MAX,CON,HL I , HL I R , HS F , F , M I  I , MPS , NNMM , M I  A,  PPS,  ALPHA) 
NL=J J-l 

KP  =  0 

SSB=( RMWP-1.0)/ ( RPWP-1.0) 

DO  302  NJ=1,NI 

SHEAR  WALL  ANALYSIS  BEGINS. 


IK  =  0 
IX  =  0 
MMMM-1 

DO  131  K=1,NS 
H( NJ,K)=F( 1 , K ) 

131  CONTINUE 

132  CALL  BAKA  ( I K , I  X , M MM M , L , NS , ND , EW , KC , KB, HS , DW , PS , NNM , MPF , MPW , 

1MPC,SC,SBW,SBF,HBW,MAX,C0N,HSF, F,MI I , NNMM , T MOM, BMOM , MOM W , ROT B , TMOM 
2F, BMOMF ,RQTF, DEFF, ROT FF, ROTO, ITER, NJ,DMCB , DMCT, DMBF , DMF B , DMB W , 

3DMW B, MMCB, MMC T, MOMBW, MOMWB, FOC, SHEC, C SHEC , SHE ARW, NCYCLE, SHEAR, COEF 
4F,CR0TF  ,ZM,MIA) 

IF(  IX  . EQ.  1 )G0  TO  245 
IF (  IK  .EQ.  1 )  GO  TO  272 

233  IF(PPS(1)  .LE.  0.0)  GO  TO  260 

I F ( MMMM  .EQ.  1  .AND.  NJ  .EQ.  1)  GO  TO  235 
I F ( MMMM  .EQ.  1)  GO  TO  256 
DO  234  N=1 , NS 

IFiABSI ICDEFF(MMMM,N)-CDEFF(MMMM-1,N) )/ CDEFF ( MMMM , N ) )  .GE.  CON  .OR 
1.  ABS( ( CROTF ( MMMM ,N ) -CROTF ( MMMM-1 , N) ) /CROTF ( MMMM, N ) )  .GE.  CON)  GO 
2T0  251 

234  CONTINUE 
GO  TO  260 


ooo  ooo  noon  oooo 
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OUTPUT  STATEMENTS  FOR  RESULTS  OF  FRAME  ANALYSIS  BY  SUBROUTINE  ROFA 
(WITHOUT  THE  EFFECT  OF  AXIAL  LOADS) 

235)  WRITE(6,240) 

240  F0RMAT(/45X, 'RESULT  WITHOUT  THE  EFFECT  OF  AXIAL  LOAD'//) 

CALL  ROFA  < ROTFF , SHEC , FOC , MMC6 , MMCT , ROT G, NS , NC YCL E , I  TER , MMMM ) 

DO  242  N= 1 , NS 

MOMFB ( N )  =  S  BF (N) *6. 0*ROTFF ( ITER, N ) 

MOMBF ( N ) =  MOMFB ( N ) 

SHEARR(N)=2.0*MOMFB(N)/HBF(N) 

WRITE (6, 241)  N , MOMB W ( N ) ,MGMWB(N) , MOMBF (N) , MOMFB (N) , SHE ARW ( N ) , SHEAR 
1 R  (  N ) 

241  FORMAT ( 1 3X , I  3 , 1 X , 4F 19 . 2 , 5 X, F 10 . 2 , 7 X , F 1 0 . 2 / ) 

242  CONTINUE 
WRITE! 6,21 ) 

DO  244  K= 1 , NS 

F 1 1  ( K ) =F ( 1 , K ) — FOC ( K  ) 

F ( 1 ,K )=F1 1 { K ) 

244  CONTINUE 
IX=1 

GO  TO  132 

245  WRITE<6,22) 

OUTPUT  STATEMENTS  FOR  RESULTS  OF  SHEAR  WALL  ANALYSIS  BY  SUBROUTINE 
SR3  (WITHOUT  THE  EFFECT  OF  AXIAL  LOADS) 

CALL  SR3  ( BMOMF , TMOMF , RGTB  »NS  »  F 1 1 , FOC  *  SHE AR , ROT  F»  DEFF  »  CROTF , CDEFF , 
1MMMM ) 

IX  =  0 

WR I T  E ( 6 , 2 1 ) 

ADDITIONAL  HORIZONTAL  LOADLTO  SIMULATE  P-OELTA  EFFECT 

251  I F { MMMM  .EQ.  MAX)  GO  TO  254 
256  PPS ( NS+1 )=0. 0 

HS(NS+1 )  =  5000 . 0 

CDEFF ( MMMM, NS+1 ) =CDEFF ( MMMM, NS ) 

F( 1 ,1)=H(NJ, 1 )+ ALPHA* (PPS ( 1)*CDEFF( MMMM, 1 ) /HS( 1 ) -PPS ( 2 ) * ( CDE FF ( MMM 
1M, 2) -CDEFF (MMMM,  1) )/HS(2)  ) 

I F ( N S  .EQ.  1)  GO  TO  253 
DO  252  N=2 , NS 

F( 1 ,N)=H( NJ,N) +ALPHA* (PPS (N) *! CDEFF (MMMM, N ) -CDEFF ( MMMM, N- 1 ) ) /HS ( N) 
1-PPS(N+1 )* (CDEFF ( MMMM,N+1 ) -CDE F F ( MMMM , N )  ) /HS (N+l  )  ) 

252  CONTINUE 

253  MMMM-MMMM+ 1 
GO  TO  132 

OUTPUT  STATEMENTS  FOR  RESULTS  OF  FRAME  ANALYSIS  BY  SUBROUTINE  ROFA 

254  wRITE(6,255) 

255  FORMAT! 10X, 'CONVERGENCE  (CALCULATION  OF  THE  EFFECT  OF  AXIAL  LOAD) 
1WAS  NOT  ENOUGH'//) 

260  CALL  ROFA  (ROTFF , SHEC ,FOC , MMCB , MMCT , ROTO , NS , NC YCL E , I TFR , MMMM ) 

DO  270  N=1 »  NS 

I F ( N J  .EQ.  1)  GO  TO  261 

IF ( A8S( DMFB( N) )  .GE.  MPF(N))  GO  TO  262 

261  MOMFB (N)=S6F(N)*6.0*R0TFF( ITER,N) 

GO  TO  263 


(  H  n:  m  k)M. 


o  r>  o  ooo  oon  noon 
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262  MOMF  B ! N  )  =  DMF  B ! N ) 

263  MOMRF(N)=MQMFB(N) 

264  SHEARR!N)=2.0*MUMFB  (N)  /HBF(N) 

WR  ITE! 6,241 ) N, M0M8W 1 N )  , MOMWB(N) , MOMBF { N ) , MUMF8 { N ) , SHE  ARW  ( N ), .SHE  ARE 

1  (N) 

270  CONTINUE 

DO  271  K= 1 , NS 

DELTA! NJ,K)=CDEFF ( MMMM , K ) 

Fil(K)=F!l,K|-FOC(K) 

F ( 1 »  K) =F1 1 ( K ) 

271  CONTINUE 
I  K  =  1 

GO  TO  132 

272  WRITE(6»21) 

OUTPUT  STATEMENTS  FOR  RESULTS  OF  SHEAR  WALL  ANALYSIS  BY 
SUBROUTINE  SR3 

CALL  SR 3  ( BMGMF , TMCMF , ROT B, NS , F 1 1 , FOC , SHE AR , ROTF, DEFF , CROTF , CDEFF , 

1MMMM) 

WRITE! 6,21 ) 

DETECTION  OF  HINGES  IN  SHEAR  WALL. 

IF1NL  .NE.  JJ)  GC  TO  273 

IF ( ASS! DELTA! NJ , NS) )  .GT.  STUPD)  GO  TO  960 

273  SSB=(RMWP-1.0)/( RPWP-1.0) 

DO  284  J  =  1 ,  L 

Dl  J)=!TMOM( JI+BMOM! J) ) /  ( 2 . 0*MP S  !  J  )  ) 

D1!J)=D!J)*MPS(J) 

IF(ABS(D(J) )  .LT.  1.0)  GO  TO  284 
I F ( NL  .NE.  JJ)  GO  TO  274 
CPHI t J)=(D! J)-1.0)/SSB  +  1.0 
Mil! J )=D ( J )  *  MIA(J)  /  CPHI ! J ) 

274  KP  =  K P+ 1 

I F ( KP  .GE.  2)  GO  TO  281 

275  WRITE(6,280) 

280  FORM  AT ( 45 X, 3 7H DETECTION  OF  HINGES  IN  THE  SHEAR  WALL///) 

281  WRITE(6,282) J,Dt J) ,01! J) 

282  FORMAT! 5X,20HHINGE  IN  SECTION  NO. , I  3, 2X , 5 1HR AT  I  0  OF  MOMENT  IN  WALL 
1  TO  PLASTIC  MOMENT  CAP AC  I TY= , F 5 . 2 , 2X , 1 5HM0M ENT  IN  W ALL= , F14. 2/ ) 

WR  I  TE(  6 ,283)  M I  I ( J ) , M I  A l J ) 

283  FORMAT  (5X, 'REDUCED  MOMENT  OF  INERTIA  = • , F 13. 2 , 3X , ' I NI T I AL  MOMENT 
1  OF  INERTIA  =• , F 1 3 , 2/ / ) 

284  CONTINUE 
WRITE! 6,21 ) 

DETECTION  OF  HINGES  ON  FRAME  BY  SUBROUTINE  SRI. 

CALL  SRI  ( DMF  B ,DMBF , DMBW , DMWB  »  DMCT , DMCB, MOMFB, MOMBF , M0M8W , MO 

1MWB, MMCT, MMCB,NS,M3 , MPF ,MPW,MPC ) 

I  K  =  0 

I E ( NL  .tQ.  JJ)  GO  TO  285 
I F ( K P  .GE.  1  .OR.  M3  .GE.  1)  GO  TO  295 

HORIZONTAL  LOAD  IS  INCREMENTED 


265  I F ( N  J  .EQ.  NI)  GO  TO  960 
WR ITEl 6,21  ) 
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WRITE! 6,291 ) 

291  FORMAT! 50X, • HOR I ZONTAL  LOAD  INCREMENTED'//) 

WR I TE ! 6 , 2  92 ) 

292  FORMAT! 10X, ' FLOOR  NO . * , 5 X , • HOR I ZONT AL  LOAD ( K) 5X,  • VERT  ICAL  LOAD! 
IK) • //) 

DO  294  K=  1 , NS 
F ! 1 , K) =H( NJ, K ) +HL I *FD ! K ) 

I F ! I VL  .LE.  0)  GO  TO  296 
PPS(K)=PPS (K) +HLI*PPD(K) 

296  WRITE!6,293)K,F( 1»K),PPS(K) 

29  3  FORMAT!  13X,I3,11X,F6.2,16X,F8.2/) 

294  CGNTINUE 
KP=0 

GO  TO  302 

HORIZONTAL  LOAD  IS  DECREMENTED 

COMPUTE  SECTION  HEIGHT  FROM  BASE  IN  INCH 


295  I F ( N J  .EQ.  NI)  GO  TO  960 

ST0PD=30. 0*ABS( DELTA ( NJ , NS) ) 

WRI TE ! 6  »  300 ) 

300  FORMAT! 1HK,15X,' HORIZONTAL  LOAD  DECREMENTED  TO  MAKE  THE  FRAME  ELAS 
IT  I C  SO  THAT  HORIZONTAL  LOAD  CAN  BE  INCREMENTED  SLOWLY'///) 

WRITE! 6,292) 

DO  301  K  =  1 »  NS 

F ( 1,K)=H{NJ,K)-HLI*FD(K)+HLIR*FD(K) 

IF! IVL  .LE.  0)  GO  TO  303 

PPS  (  K)  =  PPS(K)-HLI*PPD(K)«-HLIR*PPD(K) 

303  DMFB! K ) =0 • 0 
DMBF  t  K ) =0 . 0 
DM WB ( K ) -0 . 0 
DMB  W ( K ) =0 . 0 
DMCT ( K)-0.0 
DMCB(K)=0.0 

WR I TE ! 6 ,293 )  K , F { 1 , K ) , PPS ! K ) 

301  CONTINUE 
HL  I  =  HL I R 
NL= J  J 

K  P-0 

302  CONTINUE 

OUTPUT  THE  ROTATION  AND  DEFLECTION  OF  EACH  HORIZONTAL  LOAD 

960  DO  907  N=1 , NS 
WRITE!6,915)  N 

915  FORMAT (48X, • LOAD-DEFLECTION  DATA  FOR  FLOOR  NO.  *,I3//) 

WR I T  E ( 6 , 90  5 ) 

905  FORMAT!  10X, 'HORIZONTAL  LO AD ( K I P ) ' , 5 X , • DE FL FCT I  ON ( I N ) ' / / 1 
DO  906  K  =  1 , N  J 

WR I TE ( 6 , 908 )  H ( K , N ) , D ELT A ( K , N ) 

9C8  FORMAT! 16X,F8.2, ilX,E13.6/ ) 

906  CONTINUE 

907  CONTINUE 
400  CONTINUE 

STOP 

END 


. 
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SUBROUTINE  SR  READS  IN  DATA  -  COMPUTES  CO-ORDINATES  OF  EACH 
FLOOR  FROM  THE  BASE  AND  THE  STIFFNESS  OF  MEMBERS  -  ALSO  ALL 
THE  INPUT  QUANTITIES  ARE  PRINTED  OUT 


SUBROUTINE  SR  ( L , NS , NO ,  EW , KC , KB , FD , N I , HS , DW , P S  ,  MPF , NNM , 

1  MPW » MPC »  SC,SBW,SBF,IVL,PPD,  HBF » HBW  , 

2  RMWP,RPWP  »MAX»CON,  HLI ,  HL  I R»HSF,F,MI I , MPS , NNMM , M I  A , PP S , A LPH A ) 

DIMENSION  FD( 30) ,HS (30) ,DW( 30  ),PS(30),  SC(30), 

1  SBF  (30)  ,SBW(30)  ,  HBF(30), 

2HBW( 30) ,HSF<30) , F(99,30) , PSW(30) ,PPS( 30) ,  PPD( 30) 

REAL  MI (30 ) ,MPF( 30) , MPW (30) ,MPC (30) ,MPSW( 30 ) »MIC( 30) »MI BF ( 30 ) , 

1  MI BW( 30) ,NM( 30) »  NNM (  300 ) , M I  I { 3 00 ) , MP S (  300 ) , NNMM( 3  00) ,MI A( 300) , 
2KB, KC 

READ (5, 200) KB, KC ,EF , E rt , NS , ND , MAX , CON , HL I » NI , HL I R 
200  FORMAT ( LX,2E18.5,2F7.0,3I3,2F5.2,I3,F5.2) 

READ(5,47C)  RM1«P,RP1«P,  ALPHA,  IVL 
470  FORMAT ( 1X,2F7.2,F5»2, 13 ) 

DO  202  K= I , NS 

READ( 5,204) FD(K) ,HS(K) , HBF (K) , HBW ( K ) , DW ( K ) , M I ( K ) , M IC ( K ) , 

1MIBF(K) »  M I B W ( K ) 

2C4  FORMAT!  1 X , F 5 . 2 , 4F 7 . 2 , 4F 1 1 .2) 

202  CONTINUE 

DO  102  K= I , NS 
READ(5,101)PSIK)  ,PSW(K) 

101  FORMAT! 1X,2F8.2) 

102  CONTINUE 

DO  205  K  =  1 , N  S 

RE  AD ( 5, 500 )  MPF (K) , M PW ( K ) , MPC ( K ) , MPSWl K ) 

500  FORMAT!  17X  , 3F  14 . 2 , F2 1 . 2 ) 

205  CONTINUE 

ALL  FORCE  UNITS  ARE  IN  KIPS  AND  ALL  LENGTH  UNITS  ARE 
IN  INCHES  UNLESS  STATED  OTHERWISE. 


L=NS*ND 
NDN=0 
K= 1-ND 
CB=0 . 0 

DO  22  J=1 , NS 
NM(J)=HS(J)/FLGAT(ND) 

HSF ( J)=CB+HS( J ) 

CB=HSF ( J) 

SB  W (  J)  =  EF*MIBW( J  )  /HBW ( J ) 
SBF( J)=EF*MIBF( J ) /HBF (J ) 
SC(J)=EF*MIC(J)/HS(J) 

F ( 1  *  J ) =FD ( J ) 

PPS( J)=PS( J)*PSW( J) 

PPD( J)=PPS( J) 

K-K+ND 
NDN=NDN+ND 
DO  23  N=K, NON 
NNM(N) =NM( J) 

Mil (N)=MI( J) 


, 
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M  I  A l N ) =M I (  J) 

MPStN)  =  MPSWJJ) 

23  CONTINUE 
22  CONTINUE 

COMPUTE  SECTION  HEIGHT  FROM  BASE  IN  INCH 


ec=o.o 

DO  25  K=1,L 
NNMM(K)=BC+NNM!K) 

BC=NNMM ( K ) 

25  CONTINUE 

InR  I  T  E  (  6,324) 

324  FORMAT (  IHI  ) 

OUTPUT  STATEMENTS  FOR  DATA  OF  FRAME  AND  SHEAR  WALL. 

WR I  T  E  !  6  »  5 1  ) 

51  FORMAT ( 50Xf29HUATA  FOR  SHEAR  WALL  AND  FRAME//) 

WRITE(6,351)  NS 

351  FORMAT ( 10X»56HNUMBEP  OF  STORIES 

1  ,13/) 

WRITE! 6,352)ND 

352  FORMAT ( 10X , 56HNUMBER  OF  DIVISION  TO  BE  MADE  IN  A  STORY  = 

1  ,13/) 

WRITE! 6»355)  Ew 

355  FORMAT ( 10X, 53HM0DULUS  OF  ELASTICITY  FOR  WALL 
1  »  F 9 . 2/  ) 

WRITE! 6,356)  EF 

356  FORMAT ( 10X , 53HMODULUS  OF  ELASTICITY  FOR  FRAME 
1  »  F  9  •  2  /  ) 

WR I T  E ( 6  »  354 )  KC 

354  FORMAT! 10X, 43HSPRING  CONSTANT  AT  BASE  OF  FRAME  ^  ,E19.5/) 

WRITE!6,353)  KB 

353  FORMAT ( 10X , 43HSPR I NG  CONSTANT  AT  BASE  OF  WALL  =  ,E19.5/) 

WRITE!6,510)  CON 

510  FORMAT! 10X, 'CONVERGENCE  L I  MIT*  , 24X , « =  * , F2C . 5/ ) 

WRITE!6,501)  ALPHA 

501  FORMAT! 1  OX, • ALPHA'  , 36 X , '  =  •  , F 20 . 3/ / ) 

WRITE! 6,364) 

364  FORMAT!  1HK) 

WRITE! 6,357) 

357  FORMAT!  10X , 10HFL00R  NO . / , 2X, 10 HFORC E ( K I P ) , 3 X, 1 7HM0MENT  OF  INERTIA, 
13X, 16HST0RY  HE IGHT ! I N ) , 3X , 14HW ALL  W I DT H (  I N )  , 3 X ,  1 5H AX  I AL  LOAD(KIP), 
23X, 15HAXIAL  LOAD(KIP) ) 

WRITE! 6,358) 

358  FORMAT! 10X,9HSTQRY  NO .  ,  1 8 X , 1 3H0F  WALL  (  I N4 )  , 43X , 9H0N  COLUMN , 10X , 7H 
2QN  WALL//) 

DO  10  N=1,NS 

WRITE(6,359)N,FD!N) »MI (N) ,HS(N)  ,DW!N) , P  S l N) , PS W ( N ) 

359  F0RMAT!13X,I3,6X,F8.2,7X,F12.2,9X,F8.2,11X,F7.2,9X,F8.2,9X,F8.2/) 
10  CONTINUE 

WRITE! 6,366) 

366  FORMAT (  1 HK ) 

WRITE! 6,360) 

360  FORMAT!  10X, lOHFLOOR  NO . / , 2 X , 14HL ENGTH  OF  BE AM , 3 X , 1 4HL ENGTH  OF  BEAM 
1 ,3X, 17HM0MENT  OF  I NERT I  A , 3X , 17HM0MENT  OF  I NERT I  A , 3 X, 1 7 H MOME N T  OF  I 
2NERTIA) 

WRITE! 6,361) 


. 
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361  FORMAT { 10X,9HST0RY  NO.,4X,13HON  ROLLER ( I N ), 3X , 1 4HC0NNEC TED  W  I  T  H  »  4  X 
1*1 4HQF  COLUMN ( IN4  )  ,5X,  1  7HOF  BEAM  ON  ROLLER , 3X , 17M0F  REAM  CONNECTED 
2  ) 

WR  ITE(  6, 362  ) 

362  FORMAT ( 42 X, 8H*ALL ( I N ) , 32X , 5H ( I N4) , 10X , 14Hw I TH  WALL! IN 4) //) 

DO  71  J- 1 , N  S 

WRITE (6  ,3631 J tHBF( J ) ,H6W( J ) ,MIC ( J ) ,  MI BF ( J  )  , MIBW ( J  ) 

363  FORMAT (13X,l3,8X,F8.2t9X,F8.2,9X,F11.2,9X,Fll.2,9X,F11.2/) 

71  CUNTINUE 

WRITE! 6,411 ) 

411  FORMAT! 1HK,30X , • PLASTIC  MOMENT  CAPACITY  OF  BEAMS,  COLUMNS  AMD  WALL 
1  •  ///  ) 

WR ITE ( 6,399) 

399  FORMAT (  1CX, lOHFLOOR  NO./,2X, 

1  15HPLASTIC  MOMENT, 2X, 15HPLAST IC  MOMENT , 2X, 1 5HPL AST  I C  MO 

3MENT  ,  5X*  1 5HPLAST IC  MOMENT) 

WRITE!  6,406  I 

406  FORMAT ( 10X,9HSTokY  NG.,3X, 

1  15HCAPACITY  OF  THE,2X, 15HCAPAC  ITY  ('F  THF,2X, 15HCAPACITY  OF  T 

2HE  ,  5X »  1 5HCAPAC ITY  OF  THE) 

WRITE! 6,407) 

407  FORMAT ( 2 2 X  ,  15HBFAM 

2 ON  ROLLER, 2X, 15HWALL  SIDE  BEAM , 2X , 1 5HC0LUMN  (KIP-IN),  5X,15HWALL 
3  (KIP-IN) . ) 

WRITE(6,40B) 

4  C  8  FORMAT ( 26X  ,  flH(KIP-IN 

1),  7X,8H(K  IP- IN) //) 


CO  410  K=  1 , NS 

WRITE(6,4C9)  K,  MPF(K) ,MPW(K) , MPC ( K ) , MPS W ( K ) 

4  09  FORMAT!  13X, 13 ,  6X,  F 1 4 . 2 , 3 X , F 1 4 . 2 , 3 X , F 1 4 . 2 , 2 X , E 1 8 

1.8/) 

4 1 C  CONTINUE 

WRITE! 6,385) 

385  FORMAT!  1HK,45X,39HSTIFFNESSES  (EI/L)  OF  BEAMS  AND  COLUMNS//) 

WR ITE ( 6,365 ) 

365  FORMAT (  10X, lOHFLOOR  NO . / , 4X , 1 9HST I F FNE S S  OF  COLUMN , 4 X, 1 7H ST  I F FN FS S 
1  OF  BEAM, 4X, 1 7HSTIFFNESS  OF  BEAM) 

WR I T  E ( 6 , 367 ) 

367  FORMAT! 10X,9HSTORY  NO . , 1 0 X , 8H ( K I P- I N ) , 1 IX , 1 5HC0NNFC T ED  WITH,  5X,  17 
1HCN  ROLLER (KI P-IN)  ) 

WRITE! 6,368) 

368  FORMAT! 49 X, 1 2 HWALL i K I P- I N ) / / ) 

DO  72  J= 1 »  NS 

WRITE! 6, 309) J, SC! J) ,SBW!J),SBF(J) 

369  FORMAT (13X, 13, 9X,F15.2,7X»F15.2»6X,F15.?/) 

72  CONTINUE 

WRITE ( 6,481  ) 

481  FORMAT! 1HK, 40X , 53HM0MENT  CURVATURE  RELATIONSHIP  OF  THE  SHEAR 


1  WALL///) 

WR I T  E ( 6 , 4o  2 ) 

482  FORMAT! 10X,9HP0INT  NO . , 5X , 19HR AT  I  0  OF  MOMFNT  I N , 5X , l 8HR A T 1 0  OF  CU 
1RVATURE  ) 

WR I TE ( 6 ,483 ) 

483  F0RMAT(24X, 19HWALL  TO  THE  PL  AS T I C , 5X , 1 8H I N  WALL  TU  THE) 

WRITE! 6,484) 

484  FORMAT (24X, 19H MOMENT  CAPACITY  OF  ,  5X  ,  1  8HC UR  V A TUR E  OF  THE) 
wR I T  E ( 6  »  48  5  ) 

485  FORMAT (29X,9HTHE  WALL . ,  10X, 18HWALL  AT  YIELD  PI.///) 
wRITE(6,488) 

RETURN 

END 
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V* 

c 

C  SUBROUTINE  BAKA  COMPUTES  THE  DISTRIBUTION  OF  LATERAL  LOAD  BETWEEN 

C  THE  FRAME  AND  SHEAR  WALL.  MOMENTS  AND  DEFORMATIONS  AT  EACH  SEGMENT 

C  OF  THE  WALL  ARE  COMPUTED.  THE  CONVERGENCE  FORMULA  IS  APPLIED 

C  (  EXCEPT  THE  FIRST  CYCLE  )  FOR  DEFLECTION  AND  ROTATION.  THE 

C  DEFORMATIONS  COMPUTED  BY  THE  CONVERGENCE  FORMULA  ARE  ENFORCED  ON 

C  THE  FRAME  SYSTEM.  THE  JOINT  ROTATION  OF  THE  FRAME  ARE  COMPUTED  IN 

C  SUBROUTINE  FRAME. 

C 

c 

SUBROUTINE  BAKA  ( IK , I X , M MMM, L , NS , ND , EW , KC . KB , H S , DW , P S , N NM , MPF  ,  MPW f 
1MPC, SC  ,SBW,SBF  ,HBW ,MAX ,CON,HSF, F, MI  I  , NNMM, TMOM, BMOM , MOM W , RO T 8 , T MOM 
2F , BMOMFjROTF, DEFF, ROT FF, ROTO, I  TER • NJ * DMCB . DMCT , DMBF » DMFB* DMBW * 

3DMW  B,  MMCB, MMCT, MOMBW* MOMW B*FOC* SH EC* CSHEC *  SHEAR W * NCYCL E, SHEAR, CDEF 
4F,CR0TF ) 

DIMENSION  BMOM (  300) , BMOMFI 30) ,CSHEC(30) , CDEFF <  99 , 30 ) ,  CROTF (99,3^ 
I ) ,DEF{ 300) ,DEFF( 99, 30) , DMCB { 30) ,DMCT( 30) , DMBW ( 30) , DMWB( 30) ,DW( 30) , 
2DMFB130) , DMBF I  30 ) ,F (99,30 )  ,FOC< 30) ,HSF ( 30 ) , H S ( 30 ) , HB W ( 3 0 ) , P S ( 30 ) , 
3R0T( 300 ) , ROTF (99,30), ROS ( 30) ,ROB( 30 ) , ROTF  F ( 99,30),  ROTO ( 99) , SC ( 30 
4,SBW(30),SBF(30) ,SHFC(30) , SHE AR W ( 30 ) , SHE AR < 30 ) , TMQM ( 3 0 0 ) , TMOMF ( 30) 
5 , VDEF  (30) 

REAL  KB, KC, MOMB,MCMW( 30 )  , M I  I ( 300 ) , MPC ( 30  )  , MOMP( 30)  , MPW ( 30 ) , 

1 MPF (30 )  ,MMCB(  30 )  , MMCT ( 30 ) , MOMBW ( 30 ) , MOMWBt  30) , NNMM ( 300 ) , N N M ( 300) 

132  DO  232  M=1,MAX 

MOMENTS  AT  VARIOUS  SECTIONS  (KIP-IN). 

MOMB  =0 . 0 
NDN-0 
K=1-ND 

DO  135  N= 1 , NS 
MGM6=M0MB+F( M,N > *HSF (N) 

K=K+ND 
NQN-NQN+ND 
DO  134  J  =  K» NON 
RS=0.0 

DO  133  I =N , NS 

133  RS=RS+F(M, I)*(HSFI  I  >-NNMM(J) ) 

TMCM ( J ) =RS 

134  CONTINUE 

135  CONTINUE 

I F ( I K  .EG.  I . OR .  IX  .EG.  1)  GO  TO  140 
I F ( M  .EG.  1)  GU  TO  145 

140  DO  142  N*1,NS 
MQMB=MCM8+MGMW ( N ) 

R  S  =  0 . 0 

DO  141  I=N*NS 
R  S=R  S  +  MGMW I  I  ) 

141  CONTINUE 
MOMW(N) =RS 

142  CONTINUE 
NDN  =  0 
K=1-ND 

DO  144  1=1, NS 
K=K+ND 
NDN=NDN+ND 
DO  143  J=K , NDN 


•  A 
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TMGMl J ) =TMOM { JJ  +  MOM W (  I) 

BMOM  (  J  )-=TMOM  (  J  )  -♦-  SHE  AR  (  I  )*NNM( J) 

143  CONTINUE 

144  CONTINUE 
GO  TO  150 

145  BMCM(1)=MGMB 

I F ( L  .EC.  1)  GO  TO  150 
00  146  J  =  2,L 

146  BMCM (  J ) — T  MOM ( J- 1 ) 

DEFLECTIONS  AT  VARIOUS  SECTIONS  (IN) 

150  A  =  ITMGMl 1 )+MOMB) *NNM< 1) / 12.0*EW*MI I ( l ) ) 

RQTB=M0M6/KB 

ROT (  1)— ROTB  +  A 

DE  F  (  1  )  =  ROT  B*NNM (  1)+A*NNM(  D/2.0 
I F ( L  .EC.  1)  GO  TO  152 
DO  151  J=2,L 

B=(TMOM{ J)+BM0M( J ) ) *NNM ( J ) / ( 2 . 0 *EW*M I I ( J ) ) 

ROT ( J )  =  ROT ( J- 1 )  +  B 

OEF ( J ) =DEF ( J-l )+RCT( J-l )*NNM( J ) +6*NNM ( J ) / 2 . 0 

151  CONTINUE 

MOMFNTS,  ROTATIONS  AND  DEFLECTIONS  AT  EVERY  FLOOR  LEVEL 

152  DO  153  N=1,NS 
I=N*ND 

K=  I  +  l-ND 

TMCMF ( N )-T  MOM ( I ) 

BM0MF(N)=6MGM(K) 

ROTF ( M, N)=ROT ( I ) 

DEFF ( M , N )  =  DEF ( I ) 

153  CONTINUE 

I F (  IX  . Eti.  1 )  RETURN 
I  F (  IK  .EG.  1 )  RETURN 
I  F  { M  .EG.  1)  GO  TO  155 
DO  154  N=1 1 NS 

ROTF l  M ,  N ) -ROT  F (  1 ,  N) *ROTF ( M-1,N) /(ROTF ( M- 1 , N ) -ROTF ( M , N )  ) 

DEFF (M,N)=OEFF( 1 , N ) *D EF F ( M- 1 , N ) / ( DEF F ( M- 1 , N ) -DE FF ( M , N ) ) 

154  CONTINUE 

155  RQS(  1 )  =  DEFF ( M , 1 )/HS< 1) 

I F(NS  .EG.  1 )  GO  TO  161 
DO  160  1=2, NS 

ROS( I ) = ( DEFF ( M, I ) -DEFF ( M , 1-1) )/HS( I ) 

160  CONTINUE 

161  DC  162  N= 1 »  N  S 
VDEF(N)=-ROTF ( M, N) *DW(N)/2.0 
ROB(N) =VDEF( N)/HBW( N) 

162  CONTINUE 

COMPUTE  MOMENT,  SHEAR  AND  FORCE  ON  FRAME  BY  SUBROUTINE  FRAME 

CALL  FRAME  ( MAX, ROTFF , ROTO , SC , ROS , KC , S BW , SBF , ROB, ROTF , NS , ITER, MMMM 
1  ,  M,  DMC  B »  DMC  T,DMBW,OMWB,DMBF,DMFB,MPC,MPF,MPW,NJ,CON) 

I F ( N  J  .EG.  1 )  GO  TO  180 

IF  l  ABS  (  DMCB(  1 )  )  .GE.  MPC(D)  GO  TO  181 

180  MMCB ( 1 )=-KC*ROTO ( ITER ) 

GO  TO  182 

181  MMC6( 1 ) =DMCB ( 1 ) 
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182  IFINJ  .EQ.  1)  GO  TO  183 

I  F I ABS I  DMCT  (II)  .GE.  MPCIi))  GO  TO  184 
IF ( ABSI OMCBI 1 ) )  .GE.  MPC(l))  GO  TO  185 

183  MMCTI 1  )=SC ( 1 ) *14 ,Q*ROTFF ( I  TER , 1 ) +2 . 0*ROTG ( ITER ) -6. 0*R0S { 1 )  ) 

GO  TO  192 

184  MMCTI 1 ) =DMCT ( 1 ) 

GO  TO  192 

185  MMCTI  1 )=SC  l  1 ) *1  3.0*ROTFF I  ITER, l )-3.0*R0SI  1 )  ) +0.5*DMCB( 1 ) 

192  SHECIl )=(MMCT( 1 ) +MMCBI 1) ) /HS ID 

194  I F ( N S  .EQ.  1  I  GO  TO  210 
DO  205  K— 2 , NS 

IF  INJ  .EQ.  1)  GO  TO  195 

IF! ABSt CMCBIK) )  .GE.  MPCIK))  GO  TO  196 

IF( AbS(DMCTIK) )  .GE.  MPC(K))  GO  TO  197 

195  MMCBI K)  =  SC ( K ) *{4.0*ROTFFl  I  TER , K- 1 ) +2 . C*ROTFF ( I T ER , K ) -6 . 0*RQS IK)  ) 
GO  TO  198 

196  MMCB ( K ) -DMC6 I K ) 

GO  TO  198 

197  MMCBI K)=SC lK)*(3.0*ROTFF{ITER,K-l)-3.n*ROS(K)) +0.5*DMCT (K ) 

198  IFINJ  .EQ.  1)  GO  TO  199 

IF! ABSI DMCTIK) )  .GE.  MPCIK))  GO  TO  200 
I F ( ABS( DMCBIK) )  .GE.  MPCIK))  GO  TO  201 

199  MMCTI K)=SC (K) *(4.0*RGTFF(ITER,K)+2.0*RQTFFI I  T ER  ,  K- 1 )  - 6 . 0 *ROS  I  K  )  ) 
GO  TO  202 

200  MMCT I K ) =DMCT I K ) 

GO  TO  202 

201  MMCTI K)  =  SC(K) *I3.0*R0TFF(  I T ER , K ) -3 . 0*R0S I K )  ) +  0 . 5 *OMCB I K  ) 

202  SHECIK)=IMMCT(K)+MMCB(K) )/HSlK) 

204  F0CIK-1 )=SHfcC(K)-SHEC(K-l) 

205  CONTINUE 

210  FOC(NS)=-SHECINS) 

SHEAR! 1 )=0.0 

DO  219  K=1 , NS 

IFINJ  .EQ.  1 )  GO  TO  211 

IFIABSIDMBWIK) )  .GE.  MPWIK))  GO  TO  212 
IF! ABSI DMWBIK) )  .GE.  MPW(K))  GO  TO  213 

211  MOMBWl K)—SBWIK) *(4.0«R0TFFI I T ER , K ) +2 . 0*ROTF I M, K ) -6.0*R0R<K) ) 

GO  TO  214 

212  MOMBWl K )=DMBW I K) 

GO  TO  214 

213  M0MBW(K)  =  S6W(K)*(3.0*R0TFF( I T ER , K ) -  3 . C*ROB I K ) ) +0. 5 *0MW8 I K ) 

214  IFINJ  .Eg.  1 )  GO  TO  215 

IF  I ABSI DMWBIK)  )  .GE.  MPW(K))  GO  TO  216 
IFIABSIDMBWIK))  .GE.  MPW(K))  GO  TO  217 

215  MOMWBIK  )=  SB W I K ) *  I  2 . 0 *ROTFF 1 1 T ER , K ) +4 . 0*ROTF I M , K ) -6 . 0*R0B I K )  ) 

GO  TO  218 

216  MOMWB(K)=DMWBlK) 

GO  TO  218 

217  MOM W B I K ) =  S B W I K ) *  I  3. 0*ROTF I M, K ) -3 . 0*RU6 I K )  )+0.5*0MBW I K ) 

218  SHEAR W I K ) = I MOMB W I K) +MOMWB I K ) J/HBWIK) 
MOMW(K)=-MOMWB(K)-SHEARW(K)*DW(K)/2.0 
FI M*1,K )=-FOC IK) 

SHEAR! 1)=SHEAR{ 1 )*F(M+1,K) 

219  CONTINUE 

IFINS  .EQ.  1)  GO  TO  227 
DO  226  K=2 , NS 

SHEAR! K)=SHEAR(K-l)-F(M+l,K-l ) 

226  CONTINUE 

227  IFIM  .Eg.  1)  GO  TO  232 
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I F ( M  .EO.  MAX)  GO  TC  221 
DO  220  N=1,NS 

IF (ABS (  ( DEFF (M,N  )-D6FF( M-1,N) J/DFFF (M,N)  I  .GE.  CON  .OR.  ABS((ROTF( 
1M»N)-R0TF( M-  1 » N ) ) /ROTF (MtN) )  .GE.  CON)  GO  TO  232 

220  CONTINUE 
GO  TO  223 

221  WRITEl 6,222)  MMMM 

222  FORMAT ( 10X, 'CYCLE  NO  .  =  ' , I  3 , 7X, • CONV ER GE NC E  (CALCULATION  ON  WALL)  M 
IAS  NOT  ENOUGH'/) 

223  NCYCLE=M 
SHEAR  (  1  )-=0.0 
DO  224  N=i,NS 

CDEFF( MMMM,N)=DEFF (NCYCLE,N) 

CROTFIMMMM, N)=ROTF (NCYCLE,N) 

IFIMMMM  .EU.  1)  CSHFC(N)=SHEC(N) 

SHEAR! 1 )=SHEAR( 1  H-F( 1,N) 

224  CONTINUE 

I F ( N S  .EU.  1)  GO  TO  230 
DO  225  K  =  2  » NS 

SHEAR(K)=SHEAR(K-1)-F( 1,K-1) 

225  CONTINUE 

230  DO  231  K= 1 1  NS 

SHEAR(K)  =  SHEAR(K)-»-SHEC  (K) 

231  CONTINUE 
RETURN 

232  CONTINUE 
RETURN 
END 
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C 

C  IN  SUBROUTINE  FRAME,  THE  JOINT  ROTATIONS  OF  THE  FRAME  FOR  A 

C  SWAYED  POSITION  (ENFORCED  BY  THE  WALL  )  ARE  COMPUTED. 

C  THIS  IS  PERFORMED  BY  GAUSS-SEIDEL  ITERATION  METHOD.  IF  A  HINGE 

C  FORMS  IN  THE  STRUCTURE  THE  JOINT  ROTATION  EQUATION  IS  MODIFIED. 

C 

c 

SUBROUTINE  FRAME  (MAX, ROT FF , ROTO , SC ,R OS , KC , SBW , SBF , ROB , ROTF , NS , I T F 
1R,MMMM,M,DMCB,DMCT, DM B W , DMW B, DMBF , DMFB , MPC , MPF , MPW , N J , CON ) 

D I  MENS ICN  RCTFF (99,30) , ROTO (99 ) , SC ( 30 ) , RO S ( 30 ) , SBW ( 30 ) , SBF (  30 ) , R  U B 
1 ( 30) , ROTF ( 99, 30 ) , DMCB ( 30) , DMCT ( 30 ) , DMBW ( 30) ,DMWB( 30 ) , DMBF ( 30 ) , 
2DMFB( 30 ) , AA(30) ,BB( 30) ,CC( 30) ,DD( 30) ,EE( 30) ,FF( 30) ,GG( 30) , HH ( 30) , B 
3A( 30) 

REAL  KC , MPC ( 30 )  , MP F ( 30 ) , MPW ( 30 ) 

DO  174  1=1, MAX 
I F {  I  .GT.  1)  GO  TO  164 
DO  163  K=1 , NS 
ROTFF ( 1 ,K)=0.0 

163  CONTINUE 

164  N=  I 

I F ( N  .EQ.  1)  GO  TO  165 
N  =  N-  1 

165  SC ( N  S+  1 )  =0  •  0 

ROS ( NS+ 1 )=RQS (NS ) 

ROTFF ( N , NS  +  1 ) =RGTFF (N,NS) 

MPC(NS+1)=500.C 

DMCB (NS+1 1=0*0 

0MCT(NS+1)=0.0 

I F ( N J  .EQ.  1)  GO  TO  10 

I F ( A8S ( DMCB ( 1 ) )  .GE.  MPC(l))  GO  TO  20 

I F ( A3S ( DMCT ( 1 ) )  .GE.  MPC(l))  GO  TO  11 

10  ROTO ( I )= ( 6. 0*SC( 1 )*ROS( 1 )-2.0*SC( 1 ) *ROTFF (N, 1) ) / ( 4. 0*SC ( 1 )  +  KC ) 

GO  TO  20 

11  ROTO ( I )  =  { 3 • 0*  SC ( 1 ) *ROS ( 1 )-0.5*DMCT( 1) ) / ( 3 .0  *SC ( 1  >  +  KC) 

20  DO  170  J=1  ,NS 

I F ( N J  .EQ.  1)  GO  TO  21 

I F ( ABS ( DMBW ( J ) )  .GE .  MPW(J))  GO  TO  30 
I F ( A8S ( DMW  B ( J )  )  .GE.  MPW(J))  GO  TO  40 

21  AA{ J)=SBW( J) *(6.0*RGB( J ) -2.0*ROTF ( M, J ) ) 

EE( J)=4.0*SBW( J ) 

GO  TO  50 

30  AA ( J )=-DMBW ( J ) 

EE ( J )=0 .0 
GO  TO  50 

40  AA(J)-3.0*S6W(J)*R0B( J ) -0 . 5*DM W B ( J ) 

EE ( J)=3.0*SBW( J) 

50  I F ( N J  .EQ.  1)  GO  TO  51 

I F ( ABS ( DMBF ( J ) )  .GE.  MPF(J))  GO  TO  60 
I F ( ABS ( DMFB( J ) )  .GE.  MPF(J))  GO  TO  70 

51  BB ( J )  =  0 . 0 

FF  (  J )  =  6 . C*  SBF ( J  ) 

GO  TO  80 

60  BB ( J ) =-DMBF ( J  ) 

FF ( J )  =  0 . 0 
GO  TO  80 

70  BB( J)=-0.5*DMFB( J) 

FF( J)=3.0*SBF( J  ) 


•  '  .  .  ft  )  ) 
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80  IF ( NJ  .EQ.  I )  GG  fO  81 

IF  <  AB  S ( DMCT ( J  >  )  .GE.  MPC(J))  GO  TO  90 
IF  ( ABS(DMCB( J) )  .GE.  MPC(J))  GO  TO  100 

81  I F ( J  .EQ.  1)  GO  TO  82 

CC ( J ) =SC(J)*(6.0*RQS(J ) -2 .0*RC1TFF ( I , J- 1 ) ) 

GO  TO  83 

8  2  CC ( J ) =SC ( J )  * ( 6  •  C *ROS ( J ) -2.0*RO  TO (  I  ) ) 

83  GG ( J ) =4 . 0*  SC ( J ) 

GG  TO  110 
SO  CC  (  J  ) =-DMCT ( J ) 

GG ( J ) =0 . 0 
GO  TO  110 

100  CC( J)=3.0*SC( J) *ROS ( J )-3. 5*DMCB( J  > 

GG( J ) =3 . 0  *  SC ( J  ) 

110  I F ( N  J  .EQ.  1)  GO  TO  111 

I F I A8S ( DMCB ( J+ 1 ) )  .GE.  MPCIJ+1))  GO  TO  120 
IF ( ABS( DMCT( J+l) )  .GE.  MPC(J+1))  GO  TO  130 

111  U13(  J)=SC(  J  +  l  ) *( b.0*RGS(J+l)-2.0*ROTFF(N, J+l  1  ) 

HH( J) =4.0*SC ( J+l ) 

GO  TO  140 

120  DDI J)=-CMCB( J+l ) 

HH( J)=0.0 
GO  TO  140 

130  DU(J)=3.0*SC ( J+l )*ROS ( J+l ) -0.5* DMCT ( J+i ) 

HH( J)=3.0*SC( J+l ) 

1 4  C  BA(J)=EE(J)+FF(J ) +GG ( J  )  +HH ( J ) 

IF(AQS(8A(J) )  .LT.  0.0001)  GO  TO  170 
ROTFFCI »J)=(AA(J)+BB(J)+CC(J)+DU(J))/BA(J) 

170  CONTINUE 

171  ITER=I 

IFIITER  .EQ.  1IGC  TO  174 
I F ( A6S( OMCB( 1 ) )  .GE.  MPC(l))  GO  TO  176 

IF { A  B  S (  (ROTO( I ) -  ROTO (  I —  1 )  ) /ROTO ( I  ) )  .GT.  CON  )  GO  TO  173 
176  DO  172  K=1,NS 

IF{ ABS(6A(K) )  .LE.  C.0001)  GO  TO  172 

IFIABSUROTFFU  ,K)-ROTFF(  I-1*K)  J/ROTFFI  I,K)  )  .GT.  CON  )  GO  TO  173 

172  CUNTINUE 
RETURN 

173  IF ( I  .EQ.  MAX )  GU  TO  175 

174  CONTINUE 

175  WRITE(6,181)  M,MMMM 

181  FORMAT!  10X,  'CYCLE  NU.-= '  » I  3  »  2X ,  '  TO'  ,  I  3  »  7  X  ♦  'CONVERGENCE  (CALCULATION 
1  ON  FRAME)  WAS  NOT  ENOUGH'/) 

RETURN 

END 


« 


noo  ooooo 
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OUTPUT  STATEMENTS  FCR  FRAME  FORCES  AND  DEFORMATIONS 


SUBROUT  INE  ROF  A  ( ROTFF , SH EC ,FOC , DMC B , DMCT , R OTO , NS , NC YCL E , I  TFR  » MMM  M 

1) 

DIMENSION  ROTFF(  99, 30)  , SH EC ( 30 ) , FOC ( 30 ) , RO TO ( 30 )  ,DMCB(  30)  ,  D”CT<  30  ) 
WR  I T  E ( 6,371 ) 

371  FORMAT ( 50X,28HRESULTS  OF  FRAME  ANALYSIS///) 

WRITE(6,372)  MMMM,NCYCLE 

372  FORM AT ( 1  OX  t  'CYCLE  NO.  MMMM= •  ,  I  3 , 2 X , ' NC YC L E= ' » 1 3 / / / ) 

WR I T  E ( 6 , 200 ) 

200  FORMAT  t 10X,7HCGLUMNS// ) 

WR ITE ! 6 , 373 ) 

373  FORMAT!  10X, 10HFL00R  NO. / , 3X , 19H JO  I  NT  ROT AT  I  ON < RAD) , 3 X , 1 8H80T TOM  MO 
1 ME NT (KIN) , 3X, 16HT0P  MOMENT ( K I N ) , 3 X , 10 HSH E AR ( K I P ) ,  3X, 10HFQRCE(KIP ) 

2) 

WRITE! 6, 374) 

374  FORMAT! 10X,9HSTORY  NO.//) 

WR!TE(6,375)RGTC< ITER) 

3  7B  FORMAT! 12X ,4HB AS E , 1  OX , E 1 3 . 6/ ) 

DO  377  K  =  1  ,NS 

WRITE (  6,37  8)K,RGTFF (  I  TER  » K ) , DMC3! KJ ,DMCT( K) , SHEC (K )  , FOC ( K) 

378  FORMAT! 13X,I3,10X,E13.6,8XfF14.2,6X,F14.2»4XfF10.2,3X»F10.2/) 

377  CONTINUE 
814  WRITE!6,386) 

386  FORMAT ( 1HK ) 

WRITE(6,201 ) 

201  FORMAT ( // L0X»  5HBEAMS// ) 

WRITE(6»380) 

380  FORMAT ! 10X,9HFLOCR  NO . , 3X , 1 6HM0MENT  IN  WALL t 3X , 16HM0MENT  IN  WA 
1LL,3X,16HMCMENT  AT  COLUMN , 3X , I6HM0MENT  AT  ROLLER* 3X» I 4H SHEAR  AT 
2THE, 3X, I4HSHEAR  AT  THE) 

WRITE! 6,381 ) 

381  FORMAT ( 22X, 16HSIDE  BEAM  AT  THE , 3X , 1 6HS I DE  BEAM  AT  THE , 3 X , 1 6H END  OF 

1  THE  BEAM  ,  3X , 1 6HEND  OF  THE  BE AM , 3X , 1 4HENDS  OF  WALL , 3X, 14HENDS 

2  OF  BEAM) 

WRITE(6,382) 

382  FORMAT (22X,16HC0LUMN  END ( K- 1 N  )  , 3X , 1 6HW ALL  END  ( K- I N ) , 3 X, 16H0N  RO 
1LLER ( K —  IN) , 3X» 16 HUN  ROLL ER ! K- I N )  , 3X , 1 4HS I OE  BE  AM ( K I P ) , 3X , 1 4H0N  RO 
2LL ER ( K I P ) / / ) 

RETURN 


END 


' 


■ 
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OUTPUT  STATEMENTS  FOR 


SUBROUTINE  SR3  (BMOMF 
1CDEFF »  MMMM ) 

DIMENSION  F 1 1 ! 30  )  , FOC 
1CR0TF(99,30) ,CDEFF( 99 
WR I T  E ( 6,55 ) 

55  FORMAT! 39X,52HSHEAR  W 
1///) 

WRITE(6,312) 

312  FORMAT ( 3X, 10HFLU0R  NO 
IP ) ,3X, 15HWALL  SHEARIK 
2  3X »  14HDEFLECT  ION  (  IN)  ) 

WRITE(6,40l) 

401  FORMAT  (  3X,9HSTGRY  NO, 

WRITEI6, 20)RQTB 

20  FORMAT! 5X , 4HB A SE  , 93 X , 

DO  400  K= 1 , NS 
WRITE(6,313)  K  ,  F  1 1 ( K ) 

1DEFF ( 1 »  K ) 

313  FORMAT l 6X,  13, 7X, F8  .2, 

13X,E13.6/) 

4CC  CONTINUE 

WR I TE ( 6  ,404 ) 

404  FORMAT! 1HK , 50X , 3 1 HCHECK  ON  SLOPES  AND  DEFLECTIONS//) 

WRITE! 6,405) 

405  FORMAT! 5X,9HFLOOR  NO. ,5X, 10HSL0PE ( R AD ) ,4X, 1 4HDEFLECT ION < IN)/  /  ) 
DO  402  N—  1 , NS 

WRITE (6, 40 3)  N,CROTF { MMMM , N) , CDEFF { MMMM , N ) 

4C3  F0RMAT(8X, 13, 7X,E13.6,4X,E13.6/) 

402  CONTINUE 

RETURN 

END 

3F  COMPILATION 


SHEAR  WALL  FORCES  AND  DEFORMATIONS 

,TMOMF,ROTB, NS, Fll, FOC, SHEAR, ROTF,DEFF,CRQTF, 

(30) , SHEAR ( 30 ) ,ROTF( 99 , 30 ) ,DEFF ( 99 , 30 ) , 

,30) » BMOMF ( 30 ) ,TMOMF (30) 

ALL  ANALYSIS  AND  FINAL  SLOPES  AND  DEFLECTIONS 

•  / ,3X  » 1 5HWALL  FORCE ( K I P ) , 3X , 1 6HFR AME  FORCE(KI 
IP ) ,9X, 17HWALL  MOMENT ! K- l N ), 9X , 1 OH  SLOPE ! RAD )  , 

,62X,6H BOTTOM, 13X,3HT0P//) 

E  13  .  t>/  ) 

♦  FOC IK) , S  HE AR ( K ) , BMOMF ! K) , TMOMF ! K ) ,ROTF( 1,K) , 
11X,FB.2,9X,F10.2,5X,F14.2,3X,F14.2,4X, El  3. 6, 


■ 


o  o  ooooon 
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IN  SUBROUTINE  SRI  FORMATION  OF  PLASTIC  HINGE  IN  A  MEMBER  OF 
THE  FRAME  IS  DETECTED 


SUBROUTINE  SRI  (  CMF B , DMBF , DMBW , DM WB , DMCT , DMCB , MOMF B , MOMBF , MOMBW , MP 
1MWB,MMCT,MMCB, NS, M3 , MPF , MPW , MPC ) 

DIMENSION  OMF  B ( 3  0 ) , DMBF (30 ) ,  DMBW (30 ) ,DMWB( 30) , DMCT ( 30 ) , DMCB (30) 
REAL  MGMF  B ( 3  0 ) , MOMBF ( 30) , MOMBW (30 ) ,  MOMW  3 ( 30 ) , MMCT (  30)  , 

1MMCB(30) , MPF (30) , MP W ( 30 ) , MPC ( 30 ) 

M  3  =  0 

DO  403  K= 1 »  NS 
0MF6(K)=M0MFB(K) 

DMBF ( K ) =MGM8F (K) 

DMBW ( K ) =MOMBW ( K ) 

DMWB(K)=MOMWB(K) 

DMCT ( K ) =MMCT ( K  ) 

DMC8( K)=MMCB( K) 

403  CONTINUE 


633 


DO  610  K=1,NS 
WRITE(6,633) 

FORMAT ( 1HK) 

IF ( ( ABS(DMCB(K) ) ) . LT . MPC ( K ) ) GO  TO  616 


M3=M  3+ 1 


WRITE(6»617)K,CMCB(K) 

617  FORMAT { 10X,50HHINGE  AT  BOTTOM  POINT 
13,5X,40HMGMENT  AT  BOTTGM  POINT  OF 
I F { DMCB ( K)  -LT.  0.0)  GO  TO  10 
DMCB( K)=MPC(K) +0.0001 
GO  TO  616 

10  DMCB ( K ) =— MPC ( K)-0.0001 
616  IF ( ( ABS( DMCT( K) )  )  .LT.MPCt K) )GC  TO  602 
M3=M3+ 1 

WRITE(6,618)K,DMCT(K) 

6 1 B  FORMAT ( 10X,50HHINGE  AT  TOP  POINT 
13,  5X»40H MOMENT  AT  TOP  POINT  OF 

I F ( DMCT ( K )  .LT.  0.0)  GO  TO  20 
DMCT(K)=MPC(K)+C.0001 
GO  TO  602 

20  DMCT ( K ) =-MPC ( K)-C.0001 
602  I F ( (ABS(DMFBIK) ) J.LT.MPF(K) ) GO  TO  604 


OF  COLUMN  IN  STORY  NO. 
COLUMN  =»F14.2»5H  K-IN/) 


OF  COLUMN  IN  STORY  NO. 
COLUMN  =,F14.2,5H  K-IN/) 


I 


I 


M  3=M  3+ 1 

WRITE ( 6,603 )K,DMFB(K ) 

603  FORMAT ( 10X»50HHI NGE  AT  ROLLER  END  OF  BEAM  ON  ROLLER  OF  FLOOR  NO.,  I 
13,5X,40HMOMENT  AT  ROLLER  END  OF  BEAM  ON  ROLLER  =,F14.2,5H  K-IN/) 

I F ( DMFB ( K )  .LT.  0.0)  GO  TO  30 

OMFB(K)=MPF(K)+O.OOOl 

GO  TO  604 

30  DMFB(K)=-MPF(K)-0.0001 

604  IF( (ABS(DMBFIK) ) ).LT.MPF(K) )G0  TO  606 
M  3=M  3+ 1 

WRITE(  6,605)K,UMBF (K) 

605  FORMAT ( 10X,50HHINGE  AT  COLUMN  END  OF  BEAM  ON  ROLLER  OF  FLOOR  NO., I 
13, 5X,40HM0MENT  AT  COLUMN  END  OF  BEAM  ON  ROLLER  =,F14.2,5H  K-IN/) 

IF(DMBFIK)  .LT.  0.0)  GO  TO  40 


, 


;  •  [ 
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DMBF(K)=MPF(K)+0.0001 
GO  TO  606 

40  DM6F ! K)=-MPF ( K)-0.0001 

606  IF ( ( ABStDMBW! K) ) ) .LT.MPW(K) ) GO  TO  608 
M  3  =  M  3+ 1 

WRITE(6,607)K,DMBW!K) 

607  FORMAT ! 10X ,50HHI NGE  AT  COLUMN  END  OF  WALL  SIDE  BEAM  OF  FLOOR  NO.,  I 
13, 5Xf 40HMGMENT  AT  COLUMN  END  OF  WALL  SIDE  BEAM  =,F14.2,5H  K-IN/) 

IF ( DMBW ( K  )  .LT.  0.0)  GO  TO  50 
DMBW! K)=MPW( K) +0.0001 
GO  TO  608 

50  DMBW!K)=-MPW(  Kl-0.0001 

608  IFI (ABS(DMWBIK) ) ).LT.MPW(K))GO  TO  610 
M  3  =  M3  + 1 

WRITE! 6,609)K,DMWB{ K) 

609  FORMAT! 1GX,50HHINGE  AT  WALL  END  OF  WALL  SIDE  BEAM  OF  FLOOR  NO., I 
13, 5X,40HM0MENT  AT  WALL  END  OF  WALL  SIDE  BEAM  =,F14.2,5H  K-IN/) 

iF(DMWBlK)  .LT.  0.0)  GO  TO  60 
DMWB! K)=MPW(K)+0.0001 
GO  TO  610 

60  DM WB { K ) =— MPW ( K)-0.0001 
610  CONTINUE 
RETURN 
END 


■ 
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